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Resin Bonding and Strength Development in 
Offset Papers 


By Charles G. Weber, Merle B. Shaw, Martin J. O'Leary, 


and Joshua K. Missimer 


\ series of experimental offset papers was made in which the strength and resistance to 


surface pick were developed mainly by bonding the fibers with a synthetic resin instead of by 


the conventional beating 


commercial Vv 


The papers were made from three widely different furnishes of 


vood pulps, and the principal bonding agent was melamine-formaldehyde resin 


In all instances, the resin-bonded papers were superior with respect to curling, oil absorption, 


folding endurance, and resistance to tear to gel-bonded control papers with comparable burst- 


ng strength and resistance to pick. 


Expansion also was improved in sulfite-soda papers 


Only small amounts of resin were required to greatly improve the strength of papers composed 


principally of short-fibered filler pulps 
I. Introduction 


‘imum strength in wood-fiber papers is cus- 
v developed by beating with the formation 
like bonds formed by the hydration asso- 
with beating. For the best results in 
printing, paper must be made with a mini- 
of beating, because the formation of gel on 
libers must be carefully controlled. The 
pment of maximum strength by hydration 
ites some of the most troublesome properties 
encountered in offset printing, namely, 
\pansion, excessive curling, slow oil absorp- 
nd low opacity 
xperumental manufacture of offset papers at 
National Bureau of Standards the use of a 
resin, widely used for imparting wet 
gth to paper, has been extended to a new 
Resin bonds between the fibers are sub- 
| for the conventional gel-like bonds formed 
tion. The resin bonding made it pos- 
‘develop suitable strength with a minimum 
effect on the other properties of the 


unum strength was developed in wood-fiber 


‘by resin bonding with only a fraction of the 
quired to develop comparable strength 


alone 
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By using resin it was also possible to develop 
suitable strength when using high contents of 
filler pulps that do not lend themselves readily to 
hydration. This is important because these 
pulps contribute to low curl and expansivity and 
good oil absorption in paper. 

It was possible to obtain much higher folding 
endurance, a very desirable property, without the 
undesirable effects that accompany the building 
up of folding strength by hydration. The results 
were secured in experimental papers made in the 
Bureau’s semicommercial-scale paper mill by us- 
ing different combinations of typical commercial 
wood pulps. A series of papers was made from 
each fiber furnished, with controlled variations 
in beating, with and without synthetic-resin 
bonding. The work was carried on with the 
counsel of an advisory committee of technical 
representatives under the chairmansbip of R. F. 
Reed, director of research, Lithographic Tech- 
nical Foundation. Previous publications result- 
ing from research work at the Bureau with the 
cooperation of the advisory committee contain 
information on the manufacture and treatment of 
offset papers for optimum results in offset printing 


|1 to 10}.' 





II. Experimental Paper-making Equipment 


The paper-making equipment at the Bureau is 
semicommercial in size, and is adapted to the ex- 
perimental manufacture of papers under condi- 
tions that simulate those of typical industrial 
plants. Detailed descriptions and photographs 
of the essential equipment are contained in pre- 
vious publications [11, 12, 13]. The equipment 
used in this particular work consisted essentially 
of a 50-lb beater with copper-lined tub and manga- 
nese-bronze bars and plate; a jordan refiner with 
bars of bronze and steel alloy, a four-plate, flat 
screen; a 29-in. Fourdrinier paper-making ma- 
chine with a wire 33 ft. long, two presses, nine 15- 
in. dryers, a machine calender stack of seven rolls, 


and a reel; and a five-roll supercalender 
III. Fibrous Raw Materials 


Inasmuch as the purpose of this investigation 
was to explore the possibilities of substituting 
resin bonding for the effects of hydration in de- 
veloping strength, only three types of commercial 
wood pulps were used. The following pulps were 
chosen: 

1. Bleached sulfite, made by cooking eastern 
spruce wood in a solution of calcium bisulfite. 
This was a relatively standard quality of pulp 
such as is widely used in offset papers. It is com- 
posed of relatively long, strong fibers that hydrate 
rather readily with beating 


TABLE | Test data on pulps used in experimental papers 
z > 2 & 
Kind ot pup | 2 | 3) 8| a | z £ 38 
. 7 3 3 . 3 ~~ 
- x 5 ¥ 5 - i ae 
i Per Per Per Pe Pur 
ent t\ cat cent ent ent pil 
Sulfite “44 6 12.0 412 ; “02 0O4 


Deciduous wood 
sulfite a3 Ff 10.9 5. 5 6.0 1. 4 ] t 


Soda ~ 16.2 $1) 15.1 15 1.8 ; aS 


2. Bleached deciduous-wood sulfite, a short- 
fibered filler type of pulp inade by cooking decidu- 
ous woods in a solution of calcium bisulfite. The 
fibers are soft and resilient and impart desirable 
resilience and absorptiveness to printing papers. 
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However, they are tender and are lin 
break up rather than to hydrate during je} 

3. Soda pulp, another filler pulp mad « }y 
ing deciduous woods in a solution of cavstie 
The fibers are characteristically soft and | 
comparable to deciduous-wood sulfite, 
what stronger and less resilient. Th 
hydrate readily. 

The chemical characteristics of these » 
given in table 1. 


IV. Manufacturing Procedure 


The general procedure comprised the mak 
papers with strengths developed by co 


variations in beating, without resin and 
synthetic-resin binders. The comparativy: 
with a number of different resins we 
investigated with the more nearly convent 
beater furnish of 75-percent sulfite-25- 
soda pulps. Subsequently other furnish, 
explored in a similar manner. 

Fifty pounds of pulp was used for each 1 


run. Before the stock was furnished to th 
in each instance, the roll was raised aboy 
bedplate a definite number of turns. The pos 
of the beater roll is expressed as the num! 
turns of the handwheel above | or below 
setting, which is the point of contact betw: 
roll and the bedplate. One turn moves ¢] 
0.008 inch. The beating pattern with resp 
sequence and time intervals of roll setti 
shown in figure 1. 

The beaten stock was dropped directly 
chest from where it was pumped via the st 
through the jordan and screen to the ma 
without the use of a machine chest. The stuf! 
was of the regulating type with constant hea: 
an adjustable orifice. The jordan was us 
mixer only, with the same setting for all pap 

Melamine-formaldehyde resin was the s 
ard bonding agent. It was added in the for 
colloidal solution made by dissolving the pow 
resin in warm water acidified with hydroc! 
acid. The solution was added contin 
uniform rate to the paper-making SLOCK 
before it entered the head box of 
machine. The temperature of the 
maintained at 90° +2° F. One percent 0 


size was used in each of the papers, he 
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HOURS OF BEATING 


| Standard practice of lowering the beater roll for 
all experimental papers 

ll papers Value for roll settings are turns 

heel above +) contact (0) of roll with bedplate One turn 
O08 ir Lighter bar down at start 


not a variable. It was precipitated with 
maker’s alum, Al (SO, 
ol the pH of the stock. 


finish imparted by the small machine calender 


, Which was used to 
No filler was used. 
elatively low; therefore, the papers were 
light supercalendering to obtain a smooth- 
omparable to commercial “machine finish.’’ 


V. Sampling and Testing 


method of sampling paper from a machine 
n the paper mill was as follows: The paper 
don the reel of the paper machine during 
lest period was slit across the web, thereby 
g converted into a stack of sheets. The stack 
then cut so as to give an entire cross section 
The resultant stack 
From the first 


sheets for sampling. 
ded into seven parts. 
sheets were laid out as the first sheets of 


wks to be built up. From the second part, 
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13 sheets were likewise distributed as the second 
sheets of the packs. From each of the seven 
parts a sheet was added to each pack. The pro- 
cedure was then repeated until all the paper had 
been distributed into the 13 packs. The odd- 
1-3-5-7-9-11-13) packs were then 


assembled, one above the other, to make one 


numbered 


bundle of paper, and similarly the even-numbered 
packs to make another. From these bundles, 
sheets were taken for the test samples. 

All physical and chemical tests of the pulps and 
papers referred to herein were made in accordance 
with the Standard Methods of the Technical 
Association of the Pulp & Paper Industry,’ except 
air permeability. This property was determined 
with a Carson Precision Permeability Tester [14], 
which measures the rate of air flow through the 
paper per unit of area with a pressure difference 
of 1 gem” The pH of mill waters at the beater, 
stuff box, and head box was determined electro- 
metrically with a quinhydrone electrode. 

A Set of the papers comprising sheets from each 
machine run was printed in six colors by the offset 
process at the Coast and Geodetic Survey to get 
information on relative printing qualities. The 
papers were conditioned in accordance with the 
best recommended practice for control of moisture 
in multicolor printing |4] and printed in a press 
room with controlled humidity and temperature 


VI. Description of Experimental Papers 


Papers were made of 75-percent-sulfite-25- 
percent-soda_ pulps, 75-percent-soda 25-percent- 


sulfite 
sulfite 25-percent-standard-sulfite pulps. The pro- 


pulps, and  75-percent-deciduous-wood 
cedure was to make control papers from each 
furnish with three different degrees of beating, and 
then repeat each control run several times with 
different amounts of bonding resin. Four synthetic 
resins were tried in exploratory machine runs. 
Since appreciably better results were obtained with 
melamine-formaldehyde resin than with the others, 
that resin was selected as the standard bonding 
agent for this investigation. 

The properties of the typical papers made are 


shown in table 2 


Available from the Association headquarters at 122 East Forty-Second 
Street, New York 17, N. Y The 


each test as we is other details of testing 


methods specify the number of spec ens 
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PaBLE 2 Properties of experimental offset papers 7 
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Ru. ons Fiber furnist ‘leas "7 ~ ~ — — - 4 Curl ethane 6 
Perce h i : Degree 
‘2 None 1 A 10 2 f ‘ HY { 22 
‘ A 2 } & ix 4' 9 - 
1 4ew 1 2 LOA is 108 1) rl 4 ‘ 
" 1 4 2 IZA 1s ue l 2 , 4.2 
4 ! 2 6A 2 2 67 tt " 7.8 t 
14 N ( 5 + 2 4 7 “4 ‘ 
148u i 2 A 0) t 24 62 7 
148 i 24 2A ‘ SS ‘ 6y 72 6.7 ‘ 
la4 2 27 A 218 ii 4 ‘ 
14 ” I7A 61 H) "4 is is 8.4 17 
1458 I None ’ 2 4A 2 ti 47 4S 7. 14 
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is ! tu N 2A , 14 17 ’ y 4 > 
1 ilf 
is i ZA t é 4 14 J 
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ist i Nor ‘ ZA 42 2 2 i" ‘ 
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1488 1 None ’ 6A 78 t ' ‘ 8 { 24 
14M ; la-2 f Nom l 12 2A t 5 ! 10 2.4 
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2 ee oe beak 28 X50:500 Pickit Fe 
¢ ¢ Perce ¢ h 
1452 ON n” 72 st Nv 0. 0084 es Wrink 
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are compared graphically. 


pretation and Discussion of Data 


objective of this work was to find the 
ivantages of developing strength by 
ng over the dev elopment of comparable 
vy “hydration.” By referring to the 
75-percent-sulfite-25-percent-soda pa- 
le 2, it will be seen that 2 percent of 
sin and only 1's hours of beating, which 
to mixing only, developed bursting 


and surface pick comparable to those 


by 9's hours of beating without resin. 
me half hours of beating with 3 percent 
achieved far better strength than was 


by hydration alone, even with the most 
ating 

lvantages of resin-bonding versus gel- 
or hydration, for developing strength in 
papers are shown graphically in figures 
In figure 2, three papers of like bursting 
The bursting 


strength of paper C was developed by 9's hr of 
beating. For paper R;, the bursting strength was 
developed with 3 percent of resin and only 1‘ hr 
of beating, and for paper F&,, with 1 percent of 
resin and 5's hr of beating. 

The graph shows the important degree to which 
resin bonding in this type of fiber furnish produces 
strength, without the objectionable effects of 
hydration on expansion, oil resistance, and in- 
creased tendency to curl. The importance of 
these properties in printing papers is well known. 

The resistance of the surface to picking, as 
measured by the wax test, is an important prop- 
erty in offset papers. The surface fibers must 
resist the strong pull of tacky inks in printing. 
Figure 3, shows the relationship between the sur- 
face pick of the sulfite-soda papers and the tend- 
ency to curl. One advantage of resin bonding ts 
demonstrated rather clearly here. 

For the papers made of 75-percent-soda—25- 
percent-sulfite furnish, resin bonding was a rela- 
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Sevent y-five-percent-sulfite-25-percent-soda papers 


Resin-Bonded Papers 





DEGREE OF WET CURL 








BA DA I2A A 6A IA 20A 22A 
SURFACE PICK - WAX NUMBER 


FIGURE 4 Relationship between sur face pick and tendene y 


fo ‘ irl for >-pere 


} ent-sul fite 25 ~ pe rcent-soda 


b papers, 
trength de elo ped by beating alone and by resin bonding 


No resin; @, resin bonded 


tively greater aid to strength development than 
with higher sulfite content. The effects of resin 


bonding with this furnish are shown in figure 4. 





Beating intervals of 5‘; hours with | 
resin, and only 1‘: hours with 2 pere: 
produced papers superior in strength, 
tion, and tendency to curl to the control , 
With this 
of furnish, hydration plays a smaller 


with maximum beating time. 


when more sulfite is present, conseq 
cessive beating has less deleterious effec: 
sion. It does, however, produce excessiy 
to about the same degree as with the paper 


dency to curl accompanying excessi\ 


taining 75 percent of sulfite. 


cannot be ascribed entirely to hydratior 
“fines’’ formed in the beater doubtless ply 
important part. 
Deciduous-wood sulfite pulp hydrates 
little during beating, and the papers conta 
75 percent of this pulp developed strength 
slowly with beating. The effects of resin bor 
It is of interest that 


again drastic beating had little effect on « 


are shown in figure 5. 
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elative prope rlies of papers containing > 5-pe recent deciduous-wood sulfite of comparable bursting strenath develo, cd 


with and without bonding resin 


which is closely associated with hydration. resin. In both instances the resin-bonded paper 
ting 9's hours without bonding resin and 5% had more desirable oil resistance and curling 
rs with 3 percent of resin produced papers with characteristics, This difference is illustrated in 
parable strength as did beating 5's hours figure 6. The straight-line relationships in this 
hout resin and 1!) hours with 1 percent of graph, as compared to the curves in figure 3, are 
doubtless accounted for by the relative absence 





of hydration in beating deciduous-wood sulfite. 


80 - 
| 


Further evidence of the different effects of 





hydration on these two furnishes are indicated in 
figure 7. This shows that  bursting-strength 





expansivity relationships fall on the same straight 





line regardless of how the strength was developed, 








DEGREE OF WET CURL 


in the instance of deciduous-wood papers, which 
do not hydrate readily. The conventional 75- 





percent sulfite 25-percent-soda papers form two 

















straight lines, parallel but separated by approxi- 





8a 10a 12a 


ately S po s of burting strength. Thus a 
qunese eaneaneeen mat ly points of burting sti ngth hu 
a ; resin-bonded paper of this type with a burst of 
Effects of} wncreasing resistance to surjace put —< . . 
: , 25 points has no greater expansion than the gel- 
ndency to curl resin bonding versus hydration . : . ; 
salen bonded control paper with a bursting strength of 
percent deciduous-wood sulfite: 25-percent-standard sulfite 


esin; @, resin-bonded only 17 points. 
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% 
BURSTING STRENGTH , POINTS 


FIGURE 7 Bursting-strength expansivity relationships for 


two different types of furnishes, with and without resin 
hond ng 


@, with resin bonding 


The permeability of paper to air is an important 
property when considering printing papers, mainly 
because it determines the behavior of the sheet 
with respect to curl and oil absorption. It has 
been noted that the three types of papers dealt 
with developed curl rapidly with increased beat- 
ing, regardless of the rate of strength increase or 
apparent hydration. This is probably the result 
of “closing the sheet” or reducing the voids, 
whether it is by gel-formation, as in hydration of 
sulfite, or the formation of fines as in the instance 
of the filler pulps, soda and deciduous-wood sul- 
fite, which are not readily hydrated. The relation- 
ship between permeability to air and degree of 
curl are shown in figure 8. 

High resistance to pick is essential in offset 
printing papers. The surface fibers must resist 
the pull necessary to transfer tacky ink from the 
rubber blanket to the paper without lifting or 
being themselves transferred to the blanket. A 
fiber or bundle of fibers loosened by the pull of thre 
ink not only leaves a blemish in the printed image 


on the sheet where the “picking” occurs, but it 


clings to the blanket where it becomes moistened 


and repels ink to print white spots on succeeding 
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Ficure 9. Relation of reststance to surface pil 
strength; all experimental papers 


0, Seventy-five-percent<deciduous-wood sulfite, 25-percent-st 
no resin; J, same resin bonded 
soda, no resin; @, same resin bonded 


percent sulfite, no resin; 4, same resin bonded 


>, Seventy-five-percent sulf 


\, Seventy-five- 
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urface strength is closely allied to burst- 
th (see fig. 9), and is improved to about 
The results of 
tests are shown in the two right-hand 
f table 2. It will be seen that a pick 
ss than 6A resulted in some picking on 

and gave feeder trouble, apparently 


degree by resin bonding. 


fa “‘limpness.”’ The printing was other- 


alisfactory, except for misregister that 


nied feeder trouble. 
VIII. Summary and Conclusions 


By adding 1 to 3 percent of melamine formalde- 
resin, Maximum strength and resistance to 
picking was developed in offset papers 
a minimum of beating. The development 
gh strength by hydration or gel-formation 
is adversely such properties as expansivity, 
and oil absorption, which are particularly 
ortant in offset printing. 
The bonding resin was very effective in develop- 
sstrength in papers containing high percentages 
soda and deciduous-wood sulfite pulps, which 
not develop strength readily by beating alone. 
ttle as 2 percent of resin, with a minimum 
wating produced strength and surface bonding 
ior to any obtained without resin, even with 
maximum of beating. 
Using three different wood-fiber furnishes, con- 
rably better all-around printing papers were 
tained with resin-bonding than without, and 
indicated savings of beating time and power 
important. 
The application of resin bonding apparently 
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provides a means of using higher contents of short 
filler-types of fibers than is possible when develop- 


ing strength by beater action only. This is 


important from the standpoint of conservation of 
fibrous raw materials and availability of pulps 
The short-fibered pulp woods are not fully utilized, 
while the supply of long-fibered softwoods is 


becoming critical. 
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Metastable Transitions in Mass Spectra of Fifty-Six 
Hydrocarbons 


By Evelyn G. Bloom, Fred L. Mohler, Jonathan H. Lengel, and C. Edward Wise 


In the mass spectra of 56 saturated and unsaturated hydrocarbons, all the observable 
metastable transitions have been identified. In 362 cases, there are 32 different transitions 
involving loss of masses 2, 15, 16, 26, 27, 28, 29, 30, 42, and 44. The occurrence of a meta- 


stable transition depends on the intensity of both the initial ion peak and the final ion peak. 


I. Introduction this dissociation phenomenon to excited states of 


the initial ion, which spontaneously go to the dis- 
he mass spectrum of a hydrocarbon as ob- sociated state with a rather long “‘life’’ according 


ed in a gas analysis mass spectrometer shows to the relation 


‘sat integral mass numbers equal to the molec- 
weight of the compound and its dissociation 


lucts. In addition to these peaks there are . —— 
where n, and n, are the number in the initial state 
monly found some small and abnormally , , 
: at zero time and time ¢. Hipple * has measured 
« peaks with maxima that are in general at i 
— values of the life + by measuring relative peak 
ntegral mass numbers. Hipple, Fox, and , is oid 
: heights with different electric fields drawing ions 
ion have made a detailed experimental . “ae . . » 
from the ionization chamber and finds values of 
vy of these abnormal peaks and have given a ' . a Re 
a ' the order of 107® second. The dissociating ions 
retical explanation. They arise from ions # 
‘ ‘ are termed metastable ions, because the life of 
t dissociate after they have traversed the ion- mae 
:; . Me Siting : , the initial state is long compared to most radia- 
elerating field. If an ion of mass m, dissociates se 
je" ' tive transitions 
ediately after traversing the electric field into meat 
ss : The authors have been making systematic 
m of mass m, and a neutral particle, then the . 
ae a measurements of mass spectra of pure hydrocar- 
irent mass m, of the ion is given by the relation , : “i 
bons and tabulating results in the form of tables 


m,—m?2/m;. (1) published by the American Petroleum Institute, 
Research Project 44.4. The tables present intensi- 
use the initial ion does not dissociate at a ties of mass peaks relative to a value 100 for the 
nite point but over a range of positions, the maximum peak in the spectrum and include all 
sulting ion peak is broad and diffuse. Hipple, peaks of intensity 0.01 or greater. Spectra at 
il.' verified this theory by showing that the ionizing voltages of 50 and 70 volts are given 
etic energy of the final ion was reduced in All measurable peaks involving metastable transi- 
portion to the mass lost by measuring the tions have been recorded, and it has been possible 
pping potential applied to the ion collector that to identify the transitions involved in nearly every 
id reduce the diffuse peak to zero. Thus, if case. Mass spectra of 56 compounds have been 
on-accelerating potential is Vand the ion published to date. 
iks in half, then a stopping potential of 1/2 is The following compounds have been included in 
feient to stop the ionized fragment that has this study: All saturated isomers (40 in all) from 
wilf of the kinetic energy. They ascribed eS 
ates J. A. Hipple, Phys. Rev. 71, 504 (1947). 


ple, R. E. Fox, and E. U. Condon, Phys Rev. 69, 347 (1946 4 Catalog of mass spectral data, tables 1 to 546, American Petroleum Insti 
ple and E. U. Condon, Phys. Rev. 68, 54 (1945 tute, Research Project 44, National Bureau of Standards 
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methane through the C eights; unsaturated com- 


ethylene, propylene, four butenes, six 


pounds 
pentenes, two butadienes, and two butynes. All 


these compounds except methane show from | to 
In all, 362 have been 
recorded and the This 
paper is a survey of the detailed data that have 
been published in the API catalog (see footnote 4) 
Enough cases have been covered to permit some 


11 metastable transitions. 


transitions identified. 


general conclusions to be drawn concerning the 


occurrence of metastable transitions. 


II. Experimental Conditions 
A consolidated mass spectrometer has been 
used in this research, and spectra have been 
obtained with uniform procedures described in 
detail in a descriptive sheet published with the 
API catalog of mass spectra (see footnote 4). 


In this spectrometer, ions produced by ar 
beam at 50 or 70 volts are drawn out of t! 
tion chamber by a small field, accelerat 
ion-accelerating voltage and bent throuch 
in a homogeneous magnetic field. In o 
the spectra, the magnetic field is held cons 
the ion-accelerating voltage varied. T! 
draw-out potential is about 1 percent of 
That is, it is a 


potential varying inversely as the appar 


accelerating voltage. 


of the ion. 

The ion current is recorded by four gal 
eters with sensitivities in the ratios 1 t 
1/10 to 1/30. 
give nearly full-scale deflection for the maxim 


A compound is run at a pres 


peak on the 1/30 seale and a peak of 0.01 per 
of maximum is detectable under optimum co 
Figure 1 illustrates a part of the recor 


the n-butane spectrum showing the four galvano 


tions. 
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in decreasing order of sensitivity from 
ttom. The 43 peak is the maximum 
ere are four conspicuous metastable 
30.4, 31.9, 


and three other small ones (24.1, 35.1, 


peaks in this range (25.1, 


7.1) that require careful measurement and 
parison With other spectra for their verifica- 
The peak at 31.9 is a good example of the 
acteristic appearance of these peaks. It is 
lree times the normal peak width. The 

-at 39.2 is the largest of these, but it is incom- 

ely resolved from the large 39 peak so that the 

of the peak height is more uncertain than 
» the other cases. 

This interference between normal peaks and 
stable transition peaks occurs frequently and 
rive rise to a large uncertainty in the peak 

However, because of the abnormal width 
«metastable transition peaks a rough estimate 
eak height can often be made even when a 

normal peak Is exactly coincident with it. 
ple’s analysis (see footnote 3) shows that the 
tive height of a metastable peak will depend 
the draw-out potential, which is in this spee- 
eter a variable inversely proportional to the 
ss It also depends on the resolving power of 


spectrometer 


III. Intensity Relations 


The relation (eq 1) between apparent mass and 
nitial and final masses is theoretically suffi- 
t to determine both m,; and m,. because there 

additional condition that m, and m,; are 
molecular 


gers equal to or less than the 


ght. In practice, m, is not always known with 


ent accuracy to give a unique answer. An 
ortant guiding principle has been the fact that 
ever a metastable transition is observed m, 
correspond to relatively large Mass peaks 

ihe Ihass spectrum 
We have looked for a quantitative intensity 
The ratio R of the observed intensity at 
to the product of the intensities at m, and m,, 


R—I(m,)/I(m,) XI(m,)), (3) 


v of the order of magnitude 0.01 when J’s 
essed on the basis of the maximum peak 
to unity. One would expect variations for 
ental reasons where different values of m, 
ompared. There is, however, a nearly equal 
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range of variation when one compares the same 
transition in different compounds. The average 
value of R for all observations is about 0.01, and 
this furnishes a criterion for estimating whether or 
not observable metastable peaks are to be ex- 
pected. Metastable transition peaks very rarely 
exceed 1 percent of the maximum intensity, and 
this is consistent with the average value of R. 

The identification of transitions has been aided 
by going systematically from simple to more 
complicated molecules in each series. In the 
saturated series, C,H»,,0, it is only the parent 
peak that is uniquely characteristic of the particu- 
lar n value. The other per ks can also come from 
dissociation of heavier molecules. Correspond- 
ingly, metastable transitions keep recurring in 
successive members of the series except for cases 


where m, is the parent mass 


IV. Results 


The 362 metastable transitions observed in 55 
compounds involve 32 different transitions insofar 


If the 


transitions are classified according to the neutral 


as initial and final masses are concerned. 


mass lost, there are 10 classes involving loss of 
mass 2, 15, 16, 26, 27, 28, 29, 30, 42, and 44, or 
H., CH;, CH,, C:H2, C.H;, C:H,, C2H;, C:Hg, 
C3H,, and C3Hg. 

Loss of mass 2.—Laoss of mass 2 is by far the must 
of transition. Eight different 
transitions involving loss of mass 2 are listed in 
table 1. 
insofar as relative intensity is concerned. The 
ratio of the peak at 110 to the product of the 
peak heights at 114 and 112 is 7.3. 


frequent type 


The last metastable transition is unique 


The average 


TABLE 1 Vetastable transitions involving oss of mass 


Apparent Init Fina 


us 


Ethane, propane, butane, ethylene, 
butenes, butadienes, and butyne 

Saturated hydrocarbons Cy, to Cs 
tenes and pentenes 

Butanes and mono-olefins Cs; to ¢ 
butadienes, and 1-butyne 

Saturated hydrocarbons Cy, to Cs and 
propene, butenes, and pentenes 

Propane, butanes, 2 pentanes, hexanes 
heptanes, octanes 

Pentenes 

Ethane (30 is parent ion 


2,5-Dimethylhexane (114 is parent ion 





value for this ratio is 0.01, and values exceeding 
The peak at 112 is very small 
in all other octanes so there is nothing unusual 


0.1 are very rare 


in the failure to observe a 110 peak in other cases. 
Figure 2 shows the 110, 112, and 114 peaks in 


2,5-dimethylhexane 






































FIGURE spectrum of 2,5-dimethylherane 


howing the peak at 110 from the transition 114*-112 
Loss of mass 15.—-Seven different transitions 
involving loss of CH, have been observed (see 
table 2). 
parent molecule ion. The first transition is of 


In all these cases the initial ion is the 
interest, because in 1,3-butadiene there are no 
CH, radicals, yet the metastable transition is 
In 1,2-butadiene there is a CH, radi- 
According 


very strong 
eal but the 28.2 peak is not observed. 
to the intensity rule, one would expect to find this 
peak in both butadienes as well as in the butynes. 
Similar transitions involving dissociation of the 
parent ion are found in all the butenes and pen- 
tenes. In the saturated hydrocarbons, there are 
only scattered instances of loss of 15 by the parent 
ion and in most of these cases it occurs with ab- 
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normally high intensity. In the hept: 
octanes the nonoccurrence of this transit 
the other isomers conforms to the inten 
Failure to observe a 58.6 peak in n-hexa 
exception to this rule. Except in n-buta 
is a CH,),CH group at one or both en 
saturated molecules in all cases where th 
table transition occurs. Hipple, Fox, a 
don (see footnote 1) note that one mieht 


=o4 -_- 
‘ 


to find a transition 72*-»57+-+-15 in n 


but it is not observed. 


PasBLe 2 Vetastable transitions nvolving loss 


1,3-Butadiene 
Butenes 
Pentenes 


Butane (fig. 1 
2,3-Dimethyl butane 
2-Methylhexane 
2-Methylheptane 


NexXan 


Loss of mass 16.—-The loss of CH, from a hy 
carbon involves either a double dissociation o 
rearrangement of hydrogen atoms in the init 


ion before it dissociates. However, 35 
this have been observed in three different tra 
tions listed in table 3. The 30.4 peak of n-but: 
is the only case involving a parent ion, and 
The nhoh- 


occurrence of 29.5 in one hexane, two hepta 


occurs with abnormally high intensity 


and two octanes is consistent with the intensit 
rule, as the product of the peak heights at 


and 417 is small in these cases. 


TABLe 3. Vetastable transitions involt ng loss 0 


Apparent Initial Final 


Tass 


iss 


39 «©6 Pentenes 

41  Pentanes, 4 hexanes, 7 
16 octane 

42 a-Butane (fig. 1 


The loss of C,H, occurs 


one transition that is observed in 41 hydrox 


Loss of mass 26. 


m.=15.3 55* 297 + 26. 


This is found in all hexanes and octanes 
but one heptane and also in all butenes a 
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occurrence conforms to the intensity 


mass 27 Loss of C,H, is observed in 
ion from parent tons: 


m,.=13.5 04° —>27T + 27. 


in both butadienes and in 2-butyne. It 
norl ally weak according to the intensity rule. 
insymmetrical molecule 1-butyne does not 
a 13.5 peak, although this peak is found in 
italiene, which is also unsymmetrical and 
res a rearrangement of hydrogen atoms to 
n half 
mass 2S Loss of C,H, appears in two 
tions that recur in many saturated hydro- 
26.0 71 437+28 occurs in hexanes, 
ytanes and 11 oetanes, 
8.2 SS D7 


+28 occurs in 5 heptanes 
9 octanes. The appearance or nonappear- 
of these transitions in-general conforms to 
ntensity rule. 
f mass 29 Loss of C,H only occurs from 
ons and is observed in three transitions 
n table 4 


Vetastahble transitions involt ng loss of mass 29 


from the parent ton 


5 Pentenes 
He ptane 


Octane and 3,4-dimethylhexane 


ss of mass 30.—Loss of C,H, can only occur 
a double dissociation or by a rearrangement. 
the dissociating ion is always the parent ion, 
le dissociation is indicated. Table 5 lists the 
transitions observed. The intensity of 61.9 
ve to the product of the intensities at 114 and 


tishigh in the octanes. 


Vetastable transitions involving loss of mass 30 


Occure 
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Loss of mass 42.—Loss of C;H, occurs in two 


transitions in heptanes and octanes: 


Mq=21.7 85*-»43*+-42 in 8 heptanes and 


1l octanes: 


M,= 32.8 99° +57 42 in 17 octanes. 


The occurrence of 21.7 conforms to the intensity 
rule, but the relative intensity of 32.8 is abnormally 
high. 

Loss of mass 44.—A peak at mass 31.4 in three 
heptanes has been tentatively ascribed to the 


following transition: 
Ma=31.4 100*—>56* +- 44. 


It is observed in n-heptane, 2,4-dimethylpentane 
and 2,3-dimethylpentane. Its occurrence does not 


conform entirely to the intensity rule. 


V. Summary 


On the.basis of momentum considerations, it is 
believed that the process of ionization by electron 
collision consists in removing an electron from the 
molecule, and that the molecule ion may be left 
in a highly excited state and subsequently disso- 
ciate into various fragments. Because of the 
fact that some of the dissociation processes occur 


with a time delay of the order of 107° second, we 


have direct experimental evidence as to the 


nature of the dissociation processes. The great 
variety of metastable transitions show that mole- 
cule ions may dissociate step by step (in all cases 
where the initial ion is not the parent ion) or by 
losing a rather large mass from the parent ion. 
It is natural to assume that loss of masses CHa, 
C.H,, C.H;, and C,H, involve a single dissociation 


process, insofar as these are natural units of the 
On the other hand, loss 


hydrocarbon molecule. 
of masses CH, and C,H, necessarily involve more 
complicated processes, either double dissociation 
or a rearrangement of hydrogen atoms. 

The generality of the relation that the intensity 
of the metastable transition depends on the in- 
tensity of the initial and final ion peaks is interest- 

Consider 
Where the 


parent peaks are small, dissociation of the parent 


ing and in some respects surprising. 
transitions involving parent ions. 
ions is most probable; yet, these are the cases 
where metastable transitions are not observed, 
This is nota paradox because long life of an excited 


Many 


state implies low probability of transition. 
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transitions must take place with “half lives” 
much less than 107° second. 


Hipple (see footnote 3) has measured the half 


life of the two transitions from the parent ion in 
n-butane 


and estimated the fraction of 58* ions that are in 
the two initial states for these two transitions at 
time zero. The estimate is that roughly 9 percent 
of the 58* ions are in these two initial states. 
The 58* peak is 12.6 percent of the maximum peak 
and accounts for about 4 percent of all the ions in 
n-butane. If all are initially ions of mass 58, 
then 96 percent of the original ions dissociate with 
a very short life, 3.6 percent are stable or at least 
have a half life much longer than 10~-° second and 
0.4 percent are in two states, each of half life 
about 2x10~° second. 

The nonoccurrence of a metastable transition 
involving loss of CH, is significant, for mass spec- 
tra indicate that loss of the structural unit CH, 
must be very frequent. It is concluded that the 
transition is always short lived. 

The transitions involving loss of mass 2 are by 
far the most probable type of metastable transi- 
tions, and it is of interest that losses of mass | or 
mass 3 are never observed in hydrocarbons. Loss 


of mass 2 occurs in almost every case W 


are two large peaks separated by two 1 
with only a few exceptions. No such ¢ 
tion can be made for loss of larger mass 
It follows from the intensity rule tha 
having similar normal mass spectra 
show the same metastable peaks. The o 
of a peak at 28.2 in 1,3-butadiene and rn 
butadiene is an exception to this. Dib: 
and Mohler * have found another exceptio 
the isomers of monochloropropene. |) 
though, the metastable peaks will not be use‘ 
chemical analysis of unknown mixtures fo: 
merely reflect features of the normal spectrun 
The compilation of mass spectra and the id 
fication of metastable transitions is conti: 
and data will appear in the API catalog 
footnote 4). 


We acknowledge the assistance of Dorothy 
Thompson and Laura Williamson in making t 
of possible values of m,. V.H. Dibeler an 
Hipple have given valuable help and advy 
this project 
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‘un Determination of Small Amounts of Oxygen in 


1,120° C, 


of small amounts of oxygen in organic compounds. 


and collection of the carbon monoxide in a flask of known volume 


nM Organic Compounds 
- By William W. Walton, Francis W. McCulloch, and W. Harold Smith 


\ method based on the Unterzaucher procedure has been developed for the determination 


The method involves the pyrolysis of 


the sample in an atmosphere of helium, conversion of the oxygen compounds that are formed 


+, to carbon monoxide by passage of the products over pellets of carbon at a temperature of 


Interfering 


substances are removed by passage of the gas through a liquid-air or potassium-hydroxide 


trap before collection of the gas. 


determined through use of the NBS colorimetric indicating gel. 


oxygen can be determined readily. 


g ta I. Introduction 


vi ln investigations of natural [1]' and synthetic 
ers, Of plastics, of petroleum products, and 
r organic materials, a direct method for the 
small amounts of 
Such a method has many 


cations, for example, the determination of 


rate determination § of 


ven was needed. 


s of oxygen-containing compounds (water, 
in purified substances, or in studies of rate 
vidation. 

Living and Ligett [2] critically reviewed methods 
have been proposed to determine oxygen in 
« compounds and favor a procedure pro- 
1 


Vv Schiitze [3, 4, 5], which was later modified 
uterzaucher [6]. The essentials of the method 
he thermal decomposition of the sample in an 
sphere of oxygen-free nitrogen, passage of 
resulting products over carbon heated to 

(' to convert to carbon monoxide the com- 
ds of oxygen that have been formed, absorp- 


of interfering gases by solid potassium 
xide, and oxidation of the carbon monoxide 
bon dioxide by iodine pentoxide at a tem- 

of 118° C. An important feature of 
maucher’s apparatus is an arrangement 


by the nitrogen can be passed through the 


brackets indicate the literature references at the end of this 


ser@aeermination of Small Amounts of Oxygen 


The percentage of carbon monoxide in the collected gas is 


As little as 0.01 percent of 


furnace in reverse to flush out any air that is 
introduced with the sample. 

The accuracy claimed by Unterzaucher for his 
method has not been obtained by other workers. 
After a critical 
Aluise, Hall, Staats, and Becker {7] recommend the 
procedure as a satisfactory method for the deter- 


investigation of this method, 


mination of oxygen in organic compounds, but they 
state that the results obtained in their laboratory 
show over-all precision and accuracy somewhat 
less than those now obtained in carbon and hydro- 
gen analyses. Unpublished work by members of 
the staff of the National Bureau of Standards has 
disclosed a number of disadvantages in the use of 
iodine pentoxide in the determination of small 
is difficult to 
obtain samples of iodine pentoxide that yield low 
blanks, and frequently samples of iodine pentoxide 
that are satisfactory will suddenly develop large 
For the accurate determina- 


amounts of carbon monoxide—it 


blanks during use. 
tion of small amounts of oxygen such variability 
cannot be tolerated. 

Recently the Gas Chemistry Section of the 
National Bureau of Standards developed a colori- 
metric indicating gel for the detection and esti- 
mation of small amounts of carbon monoxide [8]. 
[t is possible by this method to detect and estimate 
less than 1 part of carbon monoxide in 1,000,000 
parts of air. This sensitivity is far greater than 
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needed for estimating the carbon monoxide formed 
in the Unterzaucher procedure even with very 
small amounts of oxygen—20 milligrams of a 
sample containing 0.01 percent of oxygen would 


vield a concentration of carbon monoxide equal 


to about | part in 300,000 if collected in a 1-liter 


flask 

A procedure was developed for the determina- 
tion of oxygen based on the Unterzaucher method 
for converting the oxygen to carbon monoxide and 
determination of the amount of carbon monoxide 
with NBS colorimetric indicating gel." 


II. Description of the Apparatus and 
Method 


A general view of the apparatus is given in 
figure 1, and the details are shown in a schemati 


drawing, figure 2 


Fiaure 1 \pparatus for the determination 


The rubber connections used by Unterzaucher 
were avoided, because water and oxygen diffuse 
through rubber Helium was substituted for 
nitrogen because of its inertness, its greater free- 
dom from oxygen, and the possibility that it would 
be more efficient in the removal of adsorbed 
oxygen 

Furnace (5, in fig. 2) is a Nichrome-wire-wound 
type, 22 em in length, and is operated at a tem- 
perature of 400° C. The tube in this furnace is 
made of 10-mm Pyrex glass tubing and is about 
$2 em long It is sealed to the U-tube through a 
two-way stopcock that allows the gases to be 
passed either through the U-tube or to the atmos- 


phere, and it is filled with small copper spirals to 
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remove any oxygen that may be in th 
Any oxide formed on the copper spirals is 
when necessary by passing a stream of | 
through the heated tube and out to the aty 
through the two-way stopcock (8, in fig. 2 

The U-tube (10, in fig. 2) has an outs 
eter of 25 mm and a height of 30 em. T! 
the U-tube. connected to stopcock x. 
with ascarites and the arm connected to 
12 is packed with anhydrous magnes 
chlorate. 

Furnace (17, in fig. 2) consists of a 
pipe, with a diameter of 20 em and a len 
em, and is closed at both ends with 
board. The transite boards are 20 mm t) 
are drilled to hold the quartz tube and 
core. This central core is an Alundum | 
in diameter with grooves of ‘\.-in. pitel 


wound with about 30 feet of 90-percent pla 


small amounts of oxryqen in organic compounds 


Ficgure 2. Schematic drawing of apparatu 
nation of small amounts of oxygen ie organic cor 
Helium tank; (2) pressure regulator }) bubble ce 

with copper spirals furnace heated to 400° ¢ 
for furnace fi jy t rmocouple for furnace 
potentiometer for Me ring the ten 
tube fille« carite and anhydro 
o-way thermocouple switch 12 
4) quartz tube packed with gas-black pellet 
hurr 17) furnace heated to rye ¢ 


former for furnace 


~) Dewar flask z 
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grou! 
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(-pert rhodium wire, 0.8 mm in diameter. 
The W ngs should be closer at the ends of the 

they are at the center. The furnace is 
packing of Santocel and_ is 
Equi- 


with a 
at a temperature of 1,120° C 
f reached, would leave only a negligible 
of the oxygen in compounds other than 
onoxide at any temperature above about 
To be sure of approaching equilibrium 
he time that the gases are in the heated 
higher temperature is needed. Too high 
temperature increases the difficulties with the 
pparatus. The temperature selected by Unter- 
weber was found to be satisfactory. 
in fig. 2) is a spiral 


The liquid-air trap (21, 
med from 5-mm glass tubing and has a diam- 


rof 4 em. and a length of 25 em. 
he collecting flask is a round flask fitted at both 
is with capillary stopcocks. The upper stop- 

k is sealed to a standard spherical ground-glass 
joint. Both 1-liter and 2-liter flasks have been 
sed. They are calibrated by filling them with 
water, draining to a mark on the lower capillary, 
and weighing the water 

The quartz tube is shown in detail in figure 3. 
It is 66.5 em long, with an outside diameter of 
Quartz to 
Pyrex graded seals are used at points 3 and 7 
The end of the quartz 


3mm and a wall thickness of 2 mm. 


shown in the drawing 
ube through which the sample is introduced is 
ground to receive a ground-glass joint (2, in fig. 3), 
which is fitted with a capillary stopcock (1, in 
fg. 3). The quartz tube is packed with carbon 
pellets (5, in fig. 3) for a length of 10 em, starting 

a point 31 cm from the capillary sidearm. 
The packing is held in place with rolls of platinum 
gauze (4 and 6, in fig. 3) or quartz wool. Carbon 
m the products of pyrolysis slowly deposits on 
platinum gauze (4, in fig. 3) and interferes 
th the free passage of helium. This deposit is 
lily removed by withdrawing the platinum 


gauze and igniting it over a burner 


x 
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Quartz tube used in the apparatus for the determi- 
nation of oryqen 
topeock: (2) ground-glass joint }) quartz to Pyrex graded seal 


gauze or quartz wor carbon pellets; (6) platinum gauze 


quartz to Pyrex graded seal 


Determination of Small Amounts of Oxygen 


Detailed directions for the use of the indicating 
gel tubes and the preparation of the necessary 
materials are given by Shepherd [8]. The in- 
dicating gel tubes are made of 7-mm Pyrex glass 
tubing and are about 12 em long. The packing 
consists of a small piece of Fiberglas tape,’ 4.5 
em of guard gel, 14 mm of indicating gel, 2.5 em 
of guard gel, and a second piece of Fiberglas tape. 
The guard gel is purified silica gel. The indicating 
gel impregnated with a 
The puri- 


gel is purified silica 
palladium and a molybdenum solution. 
fication and preparation of the materials must be 
done with great care. 

After the apparatus is assembled, the air is 
swept out of the system with a stream of helium 
and the furnaces are heated to the proper tem- 
peratures. The helium gas is then replaced with 
a stream of hydrogen to remove adsorbed oxygen 
from the carbon packing. When helium was used, 
about 1 to 2 weeks were required before a satisfac- 
tory blank was obtained, but with hydrogen a 
low blank was obtained in 24 hours. The flow 
is continued overnight and = then 


After the apparatus is 


of hydrogen 
replaced with helium. 
ready for use, a slow stream of helium is main- 
tained continuously, but the temperature of fur- 
nace (17, in fig. 2) is reduced to about 900° C 
overnight. 

The sample is weighed in a platinum boat and 
introduced through the ground joint at (13, in 
fig. 2) into the quartz tube. The stopeocks (12), 
(20), and (13) are adjusted so that a stream of 
helium flows in reverse through furnace (17) and 
out at (13). After 2 hours of flushing, stopcock 
(13) is closed and stopeocks (12) and (20) are 
changed so that the gas flows forward through the 
furnace and out through the liquid air trap. The 
Dewar flask is filled with liquid air and placed 
around the glass spiral (21). The 
flask (23) is filled with water and connected to the 
liquid-air trap. The flask (24) is lowered, and 
the lower stopcock on the collecting vessel is 
adjusted so that the rate of flow of helium is 
The sample is then care- 


collecting 


about 30 ml per minute. 
fully decomposed by slowly advancing the burner 
and the gauze mantle over the section of the 
quartz tube that contains the platinum boat. 
The ignition requires from 10 to 15 minutes and is 
repeated to insure complete decomposition of any 


s preferable if the tubes are used the same day they are packed 
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particles that may have sublimed to other parts monoxide 24 hours after the apparatus })as boo 

of the tube. The second ignition requires only assembled, but after several weeks the ank 

10 minutes. The flow of gas is maintained until reduced to about 0.0002 percent of carbon m.noxic 

the water level in the collecting flask reaches a ; 

mark on the lower capillary arm, and then both II. Results and Discussion 

the lower and upper stopcocks are closed. The Pure compounds containing less than | pereey: 

collecting flask is disconnected from the apparatus — of oxygen in the molecule are not a 

and the upper stopcock opened momentarily to Therefore, two schemes were used to c! 

allow a small amount of air to flow into the method: (1) pure compounds contain 

flask to restore atmospheric pressure—oxygen percentages of oxygen were analyzed 

introduced at this point does not interfere with mixtures of cholesterol and polystyren: which 

the determination of carbon monoxide. The — had been analyzed previously, were prepared 

temperature and pressure in the room are noted. varying ratios to yield samples of very sma 
To determine the percentage by volume of carbon oxygen content. The results are shown in tabk 

monoxide in the collected gas, the flask is con- The materials listed in table 1 were selected | 

nected through a flow meter to an NBS indicating — determine the effect of various elements on 

tube and the gas passed through the tube at a results. They include chlorine, bromine. pity. 

rate of 70 ml per minute for a definite period (10 ; 

to 180 seconds, depending upon the concentration TABLI Determination of oxygen in purified 

of CO). The color formed in the indicating tube 

is then compared visually with standard tubes pate 

prepared at the same time by passing known 

amounts of carbon monoxide through them. 

Standard mixtures of carbon monoxide and air 

are used for this purpose. Detailed directions for ene 

the preparation of standards, reagents, and tubes Do 

are given in reference [8]. From the value - 

obtained for the percentage by volume of carbon Cetyl stearate 


Do 
Do 


collected gas (which is known by previous cali- Do 
bration of the flask), the density of carbon Cietetents 
monoxide at the temperature and pressure of the Do 
room, and the weight of the sample, it is possible a 
to calculate the percentage of oxygen in the mate- 

rial. The calculation may be made with the 


monoxide in the collected gas, the volume of the 


S-(2-methoxry-5-chloro)benzylthiuronium chlor 
> ide 
following equation Do 
Do 
P t 6 o,f x Sx VXdx0.571 Do 
ercehtage Oo ; ’ 
: t.«W 


*hosphoramide, is(dibromopheny! 
Do 

where ¢, and fy are the times of flow in seconds of Do 

the standard carbon monoxide mixture and the = 

gas collected from the sample, respectively; S is stostenel palyetyoens mixture 

hw 


the percentage by volume of carbon monoxide in __, De 
the standard mixture; |!" is the volume of the Do 
collecting flask; d is the density of the carbon — 
monoxide; W’ is the weight of the sample; and the Do 
factor 0.571 is the fraction of oxygen in carbon ™ 
monoxide his compound is being prepared at the National Bureau 


Blank determinations are made by conducting for use as a micronalytical standard. Its structure is proba 
an analysis in the same manner but without a that of O-Rensyliiwentam ehiatde; however the strustan 
, — : compound is not definitely known. See Bougault and Cl 


sample. The blank is about 0.001-percent carbon rend. 218, 559 (1944). 
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and phosphorus in addition to carbon, 

ydregen, and oxygen. The values obtained in- 

these elements do not interfere with 
met! I 

Inasmuch as many 


iquid air, the liquid-air trap was re- 


laboratories do not have 


i lable 
erTcent 


wed wiil.a potassium-hydroxide trap, which re- 


ulabl 


ck the 


lic substances but not unsaturated com- 
id additional analyses were made. The 
-hydroxide trap consisted of a straight 


aves a 

wunds, & 

wtassiul 
eof Pyrex tubing, 15 mm in diameter and 18 
ong, filled with solid pellets of potassium hy- 

oxide. Table 2 shows a comparison of the 

silts obtained with this trap and the liquid-air 
». These results indicate that the potassium- 
roxide trap may be substituted for the liquid- 
trap Without loss 6f accuracy. 


Relative effectiveness of liquid-air and potassium 


hydroxide traps 


Liquid-air 


KOH trap trap 


Substance 
©% Os by % Oo by 
weight weight 
4.1 4.1 
41 4.0 
3.9 4.1 


4.0 .y 


5 -chloro) benzylthiuronium 


‘ests were made to determine the length of time 
red to sweep the system free of air after 
oduction of the sample. It was found that 
minutes of reverse flushing is sufficient. 
ever, if oxygen is adsorbed on the sample a 

time may be required as indicated by the 


If desired, 


removal of adsorbed oxygen may be accelerated 


ts for anthracene listed in table 3. 


placing the boat containing the sample near 
end of the quartz tube and slipping over this 
tion of the tube an electrically heated jacket 
(with a liquid that boils at the temperature 
ted for heating the sample during the flushing 
been made on 


determinations have 


TABLE 3. Determination of oxygen in miscellaneous 


materials 


Substance Oxygen 
by weight 
Anthracene (reverse flushing for 2 hr 0. 018 
Do O17 
Anthracene (reverse flushing for 16 hr) . 015 
Do .O12 
Diphenyl oll 
Do mie! 
Natural rubber (unextracted 
GRS rubber (unextracted) 
Butyl rubber (unextracted 
Gutta percha (unextracted 
Popcorn polystyrene 
Diesel oil (Fischer-Tropsch 
Vistanex (dried 3 hr at 90° C 
Vistanex (dried 3 hr at 90° and 
Vistanex (dried 48 hr at 150° C 


$hrat 120°C 


various materials, and a few of these are listed 
in table 3. 

The low values reported in the table for the 
oxygen content of anthracene and diphenyl were 
obtained on samples that had been carefully 
purified. Original samples of 
contained as much as 0.35 percent of oxygen, and 
it was possible to follow the purification by the 
gradual reduction in the percentage of oxygen 
The diesel oil was a sample made in 


these materials 


present. 
Germany by the Fischer-Tropsch process. 

The results obtained with Vistanex are very 
interesting. In connection with another prob- 
lem, it had been noted previously by one of the 
authors that polystyrene, when precipitated in a 
finely divided state, slowly gained in weight 
when heated at 50° C. This rate of increase in 
The results 
that its 
increases with increase in tem- 
Similar results 


weight was accelerated at 70° C. 
obtained with Vistanex clearly show 
oxygen content 
perature and period of heating. 
were obtained with dicetyl ether. 
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Infrared Emission Spectra of Flames’ 
By Earle K. Plyler and Curtis J. Humphreys 


The infrared emission spectra of a hydrogen flame and of a natural gas flame have been 
compared in the region from 1.74 to 244. Emission bands have been observed at 1.87, 2.78, 
1.374, 6.264, and 15p that are produced by molecules of H,O and CO). Another band at 
3.3164 is attributed to methane From 3.04 to 3.84 there are a number of rotational lines 
some of which may be produced by OH and others by H,O molecules. As the temperature 
of the flame was increased, the rotational lines became more intense. Five lines in the region 
from 3.4u to 3.84 agree closely in position with the predicted rotational lines of the P branch 
of the 2-1 vibrational band of OH. These lines correspond to values of K equal to 11, 12, 
13, 14. and 15. Wave numbers of some of the observed rotational lines were found to show 
a good correspondence with predicted values for K equal to 11 tc 36, for lines extending from 
9 to 24u. In addition to this series of lines, there are other lines that are produced by the 
rotational states of H.O molecule. <A discussion is included of the various methods of 


»bserving flame spectra. 


I. Introduction under different conditions of burning, a commercial 


, gas-oxygen torch was used. The gas and oxygen 
the February 1948 number of the Journal of 


Research® there is an article on the emission 
sectra from a Bunsen flame. In that work it 
sshown that nearly all the emission spectra from 
Bunsen flame were produced by H,O and CO, 


were completely mixed before burning to produce 
a smooth, narrow, quiet flame 5 in. high on a port 
about 1 mm in diameter. The emission spectra 
of a gas-oxygen flame and a hydrogen-oxygen 
flame have been measured and the results com- 


ecules. The CO, spectrum has bands in the _ 
p 1, . —_ . pared. 
ngions Of 4.44 and l5u. The CO, emission in the ; 
)u region is composed of at least three bands II. Experimental Work 
7 7 : - . (l= . s AD) " 1 ‘ 7 , ‘ 
‘ frequencies of 667, 668, and 721 cm". A Perkin-Elmer infrared spectrometer, with 


xe . . TPES . "re , . ros j ° . “Tr , ‘ , 
¢ frequencies correspond to the changes in interchangeable prisms of LiF, NaCl, and KBr, 


kwown energy leveis. The transitions “e (O1'O) . . ’ 
energy leveis. The transitions are (01 was used as the resolving instrument. Natural gas 


WO), (0270) to (O1'0), and (10°0) to (O1'O). A 
ission of the 15—-u emission of CO, molecule 
s been given in a recent report.’ Nearly all 
energy in the region of 2.7m is radiated by 


containing about 88 percent of methane was 


7 


burned in a gas-oxygen torch, and the flame was 
changed by varying the amount of oxygen. This 


HO ‘eae Bec ol caused the flame temperature to be different, with 
ra a sr gest asaepenbiampeass bog a corresponding effect on the relative intensities 
we ae eens band of water vapor at 6.26, of the bands. Some of the bands would increase 


the emission spectr ( ‘otatio ed ; —_ ; 
nission spectrum, and many rotational in intensity and others would decrease in intensity, 


s extending from 10u to 24. 
the region of 2.84 to 3.2u, there are 15 lines 
| have been attributed to the H,O molecule. 


iorder to measure the spectrum in the 3—y region 


as the temperature was increased. : 
Figure 1 shows the emission bands from the 

flame of a gas-oxygen torch. Natural gas was 

burned with QO, in considerable amounts, giving 


——— 
od , 
of this paper were reported, in part, at the meeting of the rise to a hot flame. T he zero branch of the H.O 
So ty of America in October 1947 combination band is found to be at 1.87, and 
\. Pivler, J. Research NBS 40, 113 (1948) RP1860 . ve as 
Science 107, 48 (1948 the P and R branches at 1.95 and 1.78, respec- 
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Figure 1. 


¢ spectrometer ts were 0.18 mm in width, and an LiF prism was used 


focused on spectrometer slit, 


tively are also resolved from the Q branch. How- 
ever, the resolution of the LiF prism was not 
sufficient for resolving any of the rotational struc- 
In the region from 2.84u to 3.55 there are 
When less oxygen is used 


ture. 
many rotational lines. 
in the flame the lines from 2.84 to 3.1 stand out 
more prominently, and the faint lines in the region 
from 3.lu to 3.55 disappear. The intense band 
with a sharp maximum at 3.3154 was observed 
when small amounts of oxygen were introduced 
in the burner and a large cone was present in the 
flame. The region of the cone was focused on the 
entrance slit of tae spectrometer. The zero branch 
is sharp and of good intensity. The R branch is 
observed and is of maximum intensity at 3.164, 
and is well separated from the zero branch. 

In figure 2 the emission spectra of the hydro- 
gen flame are shown for the spectral region of 1.7 
to 3.9 4. The 1.87-u4 maximum, which is the zero 
branch of a combination band of the H,O mole- 
cule, is quite intense. The P branch is relatively 
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Infrared emission spectrum of a natural gas-orygen flame. 


methane band at 3.3 ~ was measured from a gas-air flame w 


were 0.25 mm in width 


less intense in the hydrogen flame than in the gas 
flame. The greater intensity of the gas flame, 
may be caused by the radiation of the CO 
molecule in the region of 2y. 
The zero branch at 2.714 is somewhat redw 

by the absorption of the CO, and H,O molecules 
in the air. However, in figures 1 and 2 there has 
not been any correction applied for the atmos 
A 


pheric absorption. comparison of tly 





figures indicates that there is considerable emi j 
sion from CO, in this region. The R branel 
the H,O band has the greatest intensity at 2 
and the P branch at about 2.94. As shown - 
figure 2, the rotational lines it. the P branch 
the 2.74-H,O emission band are well marked i 
19 lines are observed in the region from 2.80u' : 
3.25u. The slits of the spectrometer were (1%) 
mm wide, and the flame had considerable oxy: é 
When the flame was made hotter by increasii 
the oxygen by a small amount, these lines 
less conspicuous. The flame was adjusted so 
Journal of Research @Rhrare. 








ay » lines vere of greatest contrast to the back- 
und, oad they were observed with a somewhat 
wower slit. Under these conditions it was 


nd that several of the lines had bulges on their 
was not possible to completely separate 
complex structures with the LiF prism, but 
s highly probable that many other lines could 
sbserved in this region if higher resolution were 


The curve marked 6 in figure 2 represents the 
e spectrum from 3.04 to 3.44 when the oxygen 
nereased. The greater deflections are pro- 
ed, in part, by increasing the amplication of 
thermocouple output by about 50 percent. 
der the different flame conditions the rotational 
ponents show a considerable variation in in- 
sity. At wavelengths beyond 3.1 the lines 
e greater intensity, and it has been possible to 


nd the observation of rotational structure to 


amplication by a factor of two compared to the 
value used in the region from 1.74 to 3.4u. The 
wavelengths of the rotational lines are listed in 
section IIT. 

The infrared emission in the region from 3.54 to 
5.0u from the flame was studied, when mixtures 
of natural gas and oxygen were used. It was found 
that the emitted energy could be more than 
doubled by increasing the oxygen. The correc- 
tions for the atmospheric absorption and also the 
self absorption of the flame were calculated from 
observed data. By passing radiation of one flame 
through another flame, it was found that the 
flame absorption amounted to about 20 percent 
at 4.22y. 
flame focused on the entrance slit and placing 


This was measured by having one 


another flame in the path between the first flame 
and the entrance slit. The sum of the energies 


of the two flames burning separately was compared 

















4 su. Curve ¢ in figure 2 was measured by in- with the energy of the two flames burning at the 
sing the slits to 0.25 mm and increasing the same time. The corrected energy curve for the 
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Infrared emission spectrum of a hydroge n-oxrygen flame h), 


Hotter flame and amplification increased 50 


. hotter flame and amplificat on increased 100 perce nt, 
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tu to Sw region did not vary by a significant 


amount from the curve shown in the previous 
paper (see footnote 2), and it is not shown here. 
In figure 3 is shown the 6.264 H,O emission 
band 
with a considerable amount of oxygen. 


This was observed in the hydrogen flame 
A rock- 
salt prism was used in the spectrometer, and the 
slits were 0.280 mm wide. There is very little 
emission in the center of the band, and the R 
branch has maximum energy at 5.35y, and the P 
branch at 6.754. The rotational structure is not 
well resolved, but three strong lines can be seen 
between 7.54 and 8.0u. The energy radiated by 
the flame is of low intensity. 
the record of the energy emitted from 10u to 1l5u 
The position of the 


In figure 4 is shown 


where it is a maximum. 
maximum of intensity within this region varies 
with the temperature of the flame. The energy 
emitted by the flame is very much dependent on 
the amount of oxygen that is introduced. For 
example, if the flame is hot, with a relatively 
high proportion of oxygen, the energy is about 
constant in the 10u to liu region. For cooler 
increases 


marked 


flames with less oxygen, the energy 
continuously from 104 to 15a. The 


changes in relative intensity of the rotational lines 
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URE 3 Em ss7on spect am of the @ JO hand of the 


20 molecule as obtained from a hydrogen-orygen flame 


A rock salt prism w used, and the slits were 0.28 mm in width 
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of the water vapor molecule that are 
by flame changes may prove to be a 1 
flame temperatures. Hoy 
quantitative measurements were carried 


measuring 


relation between temperature and emitt: 
All observations were made in the port 
flame above the cone, except W hen otherw 
The rapid decrease in intensity of the 
lines of the H.O molecule in the region f; 
l5u is caused largely by the absorpti 
NaCl prism and by CO, in the atmosph 
results shown in figure 4 are for a flame e 
maximum intensity at about 12. T! 
was hotter than one showing constant en 
the considerable drop in the observed 
the 14 to 15 region is caused by absorptior 
the NaCl prism. 
also observed with a KBr prism, and the envel 
form with 


This region of the spectrum y 
of the emission lines changed 
temperature. For the hottest flames, the 1 
mum intensity was observed at about 12, 
with cooler flames, the maximum intensity 


at about l5u 
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Figure 4 Rotational lines of the hyudroge n-0. 


from 9 to Lip 


5mm wide, and an NaCl prism 


The region from lou to 24u has been meas 
with a KBr prism and is shown in figure 5. Ast! 


region of 24y is approached, the energy bee 


very small. The spectrum shown in 
represents the observations on the 
flame. Its emitted energy is greater 
served for the natural-gas flame, exc: 
region of liu. 
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FIGURE 5 


Che slits were ( 


III. Discussion 


the previous paper (see footnote 2), measure- 
ts of the Bunsen flame were reported, and the 
s produced by the CO, and H,O molecules 
noted. Because of some overlapping of the 
spectra, certain regions of the spectrum could 
« well segregated. For this reason a com- 
son of the hydrogen-oxygen torch flame with 
Bunsen natural-gas—air flame would assist in 
study of the emission spectra of the H,O mole- 
This is especially true of the region 1.74 to 
The long wavelength branch of the H,O 
at 2.74 has been partially resolved, and 
it 19 lines can be seen in figure 2. By increas- 
rthe sensitivity of the instrument and increasing 
temperature of the flame, lines are observed to 
it 3.84. The lines from 3.3y to 3.8u are 
iably produced by OH and H,O molecules, 
t has been found possible only in a limited 
vion to make assignments from calculated values 
iransitions between the energy levels of either 
HO or OH. 


issed in another part of this paper. 


One exceptional instance will be 
Higher 
on will be necessary in this region in order 
t accurate and complete assignment of 


eries, 
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Rotational lines of the hydrogen-oxryqen flame from 15 to 244 measur 


Oo" : Bi 


19 20 22 


IN MICRONS 
ed with a KBr prism 


mm 2 to 24 


It has been found that these lines check very 
1 and AA=0 for 
the H.O molecule and give rise to a spacing of 


closely with the series of A} 
about 20 em However, there are several othe 
series to be expected corresponding to the higher 
rotational states of the H.O molecule, which would 
also give rise to lines in this region. So this one 
series may only, in part, explain the presence of 
these lines. The curves 6 and ¢ of figure 2 show 
the effect of temperature on the lines, and the 
relative intensities are considerably altered. In 
table 1 the wavelengths and frequencies of the 
observed rotational lines are listed. Because of 
the high temperature of the flame, the rotational- 
vibrational bands have the energy reduced near the 
center and increased at larger rotational states. 
This causes the spread in frequency from the P 
In the R branch there 
is considerable crowding of the rotational lines. 


to R maxima to increase 


The band at 3.316 yw in figure 1 was observed at 


greatest intensity when the axis of a large cone of 


the flame was focused on the entrance slit of the 
spectrometer. The R branch is readily recog- 
nized, but the P 


The zero branch is intense, and the region of maxi- 


branch of the band is weak. 


mum intensity is removed from the methane band, 
observed in absorption, by about 4 em™'. As 
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Observed rotational components in the region 


from 2.8 to 3.74 


Wave num Wavelensth ave num 


ength 
ber ber 


higher temperatures are reached in the methane 
gas there will be a gradual shift to longer wave- 
lengths for the zero branch of the band. 

The authors are indebted to G. H. Dieke of 
Johns Hopkins Univ., who furnished them with 
a table of the calculated wave numbers of the 
vibrational-rotational lines of OH and also the 
numbers for the rotational lines.* 


wave pure 


Many lines arise from the 1—>0 transition in vibra- 
tional levels and also from the 2—1 transition. 


The wave numbers calculated by Dieke were 


plotted, and it was found that the overlapping of 


the different series in the region from about 3,000 
to 4,000 em~' was so great that the available 
instruments could not separate the different states 
However, in the region from 2,500 to 3,000 em™! 
the overlapping is not so great, and a series of 
five lines that were observed in the hvdrogen 
flame check will with predicted transition in the 
2->] 

» 


shown in figure 2, ¢ 


vibrational band. The observed lines are 
With low oxygen concen- 
trations in the mixture, these lines are very weak 
The flame mixture was adjusted to give the 
greatest energy in this region, and it was a hot 


flame When 


region of 3u, it was found that the lines as shown 


observations were made in the 
in figure 2 were not as well resolved with a hot 
flame as with a moderate flame. But as the 
longer wavelength region of the band is reached, 
high flame temperatures are necessary to bring 
The 2—1 band has P, 
there are two. series 
Thus for each value of K 
with 


out the rotational lines 
Y, and R branches and 
Each line is double 
levels 


rotational there exist four lines 
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nearly the same wave number. It was no. possi) 
with the prism instrument to separate (he | 

and only one maximum could be obs 

table 2 the calculated values are ta 
Dieke’s table, and they are averaged va 


four components. The experimental 


correct to 4 em~', 


TABLE 2. Calculated and observed wave num! 
tional lines in the P branch of the 2-1 han 


P 
alculated »bserved 


These five lines that appear relatively free f; 
interference or overlapping show such close agre 
ment with predicted positions, that it app 
practically certain that they constitute th 
of the rotational structure of the P branch of 1 
2-1 band of OH, for A values, 11 to 15 inelus 
This structure is observed only because the e1 
sion of the H.O molecules has become relatiy 
weak in this limited interval. For the reg 
of shorter wavelengths, the H,O emission is 
strong as to completely mask the OH spect 
structure so that it is not possible to obser 
lines showing the characteristic intervals fo1 
2->1 band and the 1-0 band. 

In table 3 the wavelengths and wave number 
of the rotational lines from the hydrogen fla 
are listed. The observations to l5u were m: 
with an NaCl prism and from 1l5y to 244 wit 
K Br prism. 

This table comprises all of the measured va 
of the emission maxima that appeared in 
recordings of the spectrum of the hydrogen fla 
shown in figures 4 and 5. This observed spect 
appears to consist of the superposed systems 
pure rotational lines of OH and H,0. 

In order to examine the evidence that ¢! 
tional OH emission spectrum has been pro 
table 4, giving a comparison of the calcul 
observed line positions, has been construi 
calculated wave numbers, correspond! 
values 11 to 22, are obtained from Die! 


‘GG. H. Dieke, Table of OH bands, private communicat 
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hbserved rotational lines in the spectrum of the 





e ln ogen flame in the region from 9 to 24 p 
ed 
cm~ A (yp em A (yp cm 
ne tro ob- ob (ob- (ob- (ob 
served served served served served 
$0 
ues 
1.060 13. 43 745 2.11 497 
1, 038 13. 53 739 2. 41 490) 
’ 1,011 13. 70 730 20. OS 477 
G43, 13. 83 723 21.47 40) 
970 13. 99 74 21.91 156 
OH 
4 OAS 14.17 706 22. 26 449 
73 931 14. 35 697 22. 63 442 
”) 917 14. 49 ago 23. 08 453 
» 911 14. 68 6al 233.5 426) 
4 G05 14. 91 671 23. 96 417 
) RO4 15.16 659 
‘ S67 15. 35 652 
S51 | 15.61 | 641 
® a44 15. 70 637 
” S27 16. O1 62. 
2? Alf 16. 5S LS) 
“ * 41 si 16. 82 5S 
nl is sol 17.31 78 
wre ; 71 | 17.70 
PI 7Y 781 Is. 13 
t 
1@ pi 4 772 18. 54 539 
of 03 “67 18. 79 532 
M4 71 | 19.2250 
lusive 2s 753 || 19.43 1 
P eI 
alive , , : > 
\tually there are four possible lines for each K, 
e . . . . . 
nasmuch as the resolution available in these 
eriments did not permit observation of these 
et 
’ irate components, the four values are averaged 
mse! 
for ach instance to give the computed line position. 
iO } ° 
ieke’s calculations were carried out only as far 
sal » K-value, 22. We have extended them to A, 
EELLEM . . 
fla by extrapolation. It is estimated that the 
“ rin this extrapolation cannot exceed 1 em, 
h is well within the estimated experimental 
wit! 
r of the observations for this region of the 
- trum, namely 4 em 
' Of the 26 entries in table 4, all but 2 show agree- 
flay nt between the observed and calculated posi- 
‘i s within the limits of experimental error. In 
oi ost all instances of failure to find experimental 
ence for an emission line of OH, there is an 
use maximum probably originating in H,O 
F lapping the caleulated position. Some of the 
F ted lines appear only as bulges on the sides of 
1 Th profiles of other intense maxima, but agree- 
(with the predicted position seems to confirm 
rreality. The evidence for the reality of this 
‘Ystem of rotational lines of OH rests mainly on its 
pleteness rather than on the precision of the 
secrCS Girared Emission Spectra of Flames 


TABLE 


observations. 


i. 


Observed 


and 


calculated 


wave 


numbers 


pure rotational lines of OH 


calculated 


GOR 
929. 5 


95). 0 


Ou. 5 


1007.0 
1024. 5 
1041.2 


1057. 2 


> 


observed 


970 
yas 
1,011 


1, 038 
1, 060 


R 


observed 


— = 85 8S OO 
ae BEB 
<3 = 
2-38 


13. 53 
13. 03 
12. 63 
12. 3 


> 
11.85 


9. 89 


of 


the 


The appearance of a small number 


of scattered combinations dispersed through another 


complex spectrum might 
fortuitous, 


but 


when 


the 


observed 


well be regarded 
system 


as 


Is 


practically complete, the probability that it is real 


Sane 
seems to be overwhelming. 


The numerous remaining lines shown in figures 


4 and 5 and listed in table 3, are produced by higher 


rotational states of the H,O molecule. 


As men- 


tioned above. some of these overlap or coincide 
with OH, making it practically impossible to esti- 
mate the relative intensities of the lines in the 


OH system. 
intense Maximum at 
receives by 
H.O, masking an OH line in the same location 


far 


the 


A noteworthy instance is the very 


16.04, Ww hich undoubtedly 


greater 


IV. Conclusion 


contribution 


from 


The emission spectra of a gas-oxygen and hydro- 


gen-oxygen flame have been compared in the region 
Nearly all of the energy in the 


from 1.74 to 3.8u 


1.9 


region 


is produced by 


the H,O molecule. 


A small increase in the energy of the gas flame in 
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the 2-u region may be caused by emission of COs. 
The rotational lines that are observed at about 
3u are present in both flames, but they extend to 
longer wavelengths in hydrogen-oxygen flames 
The intensity of the observed bands can be in- 
creased at least twofold by increasing the oxygen. 
The emission in the region of 10u can be greatly 
flames. The 


increased with high-temperature 


energy remains about constant from 94 to l6u 
With a 


moderate flame temperature the energy increases 


and then decreases gradually to 24u 


from 9u to L5u and then decreases as the wavelength 
is increased 

In this study of the emission of flames, it has 
been shown that the amount of oxygen that is 
introduced into the flame has marked influence on 
the intensity of the band. In the 3-u region, the 
relative intensities of the rotational lines may be 
varied considerably by changing the temperature 
of the flame 


sion spectra originating in a flame, it would be 


For a quantitative study of emis- 


necessary to measure the quantity of fuel gas and 
xvgen that are mixed for burning. <A gas-oxygen 
torch or a regulated burner, giving rise to a steady 
flame, makes it much easier to obtain smooth 
traces on the recorder for the energy emitted. The 
results can be repeated closely. For quantitative 
determinations, the steady flame is very desirable 
However, in many industrial applications, the 
turbulent flame is the more common type, and 
the Bunsen burner flame approaches nearer to 


industrial burner flames. This study has been 


chiefly concerned with the study of thy 
spectra of CO, and H,O in flames. 1 
resulting from burning various hydroca: 


alcohols will also have considerable enere 
by CO, and H,O molecules. 
bands in the observed flame spectra o 


There may 


fuels that are characteristic of the mo! 
atomic groups. The bands are usual! 
intensity. 

In addition to a controlled burne: 
trometer of high aperture ratio is need 
makes it possible to observe low-intensi! 
Because of the decrease in the observed 
produced by the absorption of H,O 
molecules in the air, an instrument tha 
from atmospheric absorption is bette 
quantitative measurements. For the 
many problems, an instrument of mod 
lution is satisfactory. In the measui 
the rotational lines within a band, a gra 
trometer is better suited. 


The authors express their thanks t 
Dieke of Johns Hopkins University for suy 


calculated values of the wave numbers of th 


bands and of W. S. Benedict of the 
Bureau of Standards for advice in the int 
tion of the experimental data 
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Pressure-Volume-Temperature Data for Oxygen 
By Cyril H. Meyers 


To satisfy the demand for data on oxygen at higher pressures, an equation has been 
developed that represents the available experimental data within the accuracy of the ob- 
Values of the virial coefficients are given for densities in Amagat units in 
200° to +200° C. 


served values. 


the temperature range - Values of pressure in pounds per square inch ? 


are tabulated corresponding to temperatures from —50° to +150° F and densities up to 


320 times the density at 70° F and 1 atmosphere. These data are also presented in a chart, 


with pV/RT and pressure as coordinates. 
I. Equation of State 


use of progressively higher pressures in 
vi-pressure containers has developed a need 
lata on oxygen at pressures and densities 
Although 
quation for representing the properties of 
wen has been published by Millar and Sullivan 


rthan those at present tabulated. 


that equation ts not adequate for the present 
ose. Hence an equation has been devised 


represents all available experimental PVT 


at densities up to the critical density within 
The equation 

Exponential 
ons of temperature replace the functions 


mit of experimental accuracy. 
the Kamerlingh Onnes type. 


by that author. This equation has the form 


pV/RT B/V+C/V?2+D/V‘4, (1) 


the virial coefficients have the following 


exp M/T))- 
exp M/T)*)+C, 
T)(1—exp (—M/T)')+D, 


mission of the term in 1/V% is in accordance 
the practice of Kamerlingh Onnes. The 
titles By, By, ©), Cy, Dy, Ds, and M are con- 
*s The exponent of (17/7) is in each case, 
to (2n—1) 


practice of Kammerlingh Onnes has been 


where n is the exponent of V. 


wkets indicate the literature references at the end of this 


Data for Oxygen 


s t 4s 4 


to add terms containing higher powers of V for 
making the equation applicable at higher densi- 
ties. No attempt has been made to see whether 
the form of temperature function used here could 
be extended to such additional terms. 

A study of charts of pV versus 1/V for several 
substances shows that a linear relation exists 
between the minimum values of pV for the various 
isotherms and their corresponding values of 
density. If this linear relation may be extrapo- 
lated to low temperatures where the minimum 
occurs in the region of the subcooled liquid and 
cannot be reproduced experimentally, then since 
the density of the condensed phase is finite, the 
minimum value of pV should also be finite, even 
though negative at the lower temperatures. The 
form of temperature functions used in eq | have 
been designed such that their product with 7 
remains finite as 7 approaches zero; and in this 
respect: are considered to be better than those 
used for butadiene [2), although the constants 
for the second virial can be chosen so that over a 
wide temperature range the two functions will 
give values very closely in agreement. The 
form 2, 7*+ Bs, 


virial in the equation for butadiene, becomes at 


which was used for the third 


high temperatures a close approximation to the 
form used here. 

The value of 7 was determined by a method of 
successive approximations. The procedure con- 
sisted in assuming a value for \/, solving for B, 
and 2B, from values of the second virial derived 
from the experimental data at two temperatures, 


457 





those chosen being 0° and 100° C, and calculating 
the value of the second virial at a third tempera- 
Successive values of M7 were chosen until 


value of the second virial at the 


ture. 
the calculated 
third temperature, in this case the critical tem- 
perature, agreed with that derived from experi- 
mental data. The value obtained for M was near 
1.257. and hence the number 192.838 was adopted 
for M/. This is 154.27°K, the 
number reported by the Leiden laboratory for the 
In consideration of the law 
of corresponding 1.257, 
should be sufficiently approximate to make a good 


five-fourths of 


critical temperature. 
states, the quantity 
initial choice for fitting data to other substances, 
although the exact value five-fourths is no longer 
considered significant 

Although coefficients were determined for the 
three terms to represent all the observed data for 
oxygen within the limit of experimental accuracy, 
yet at the lowest temperatures the choice of a 
smaller value of \f would improve slightly the 
agreement between calculated and observed values 
But this change would be at the expense of exces- 
sive deviations at higher temperatures. A similar 
situation was found in attempts to apply the 
equation to nitrogen. On the other hand, the 
‘ indicate a larger 
that, 


data for nitrogen above 100° ¢ 
value of \/. This emphasizes the fact 
although eq 1 covers the entire range of data on 


oxygen, it is only empirical and is applicable to a 


limited temperature range. 
For convenience in computation, the Naperian 


number e¢ is replaced by the base 10 viving 


1—10°")+ B 


10 
1074) +D, 


numerical values to be used are as follows 


B 1.04472 B,=0.0009878 


( 647.93 10° C',= 1.426 10" 


1133.9 10 D,=5.55X 10 


273.16 
0.008223, Py= 0.001208, 


to be 125.06°h, 
ition of the data in Internationa 


und 


tained, but at 400° C the second virial 
sbout 20 percent If the exponent 
experimental data above wor ¢ rt 


t lower temperatures, 


b=0.434294(M/T) =83.7484/T 


0.434294(M/T)'=3.1143™ 10) 


d =0.434294(.M1/T )? = 43.06 > 


where 


0.434294 = loge 

and where in the right-hand side of thy 
the temperature 7 is expressed in degri 
and the volume V, in Amagat units, + 
ratio of the volume at the specified st 


and | atm. Since t! 


volume at 0° C 
pV/RT is dimensioniess in form, it 
independent of the units, and it may bi 


in any units so long as they are consist: 


II. Representation of the Experiment 
Data 


Although this work was undertaken prima 


data for use in the compress. 


calculate 
industry in the temperature range -— 50 
a much wider range has been covered 
Equation | with the constants given in t 
section of this paper represents the availabl 
within the experimental error from 14. 
200° C, at densities up to 300 or even more Ai 
units. The accuracy with which eq 1 rep 
the data is shown by figure 1, in which the co 


(1—pV/RT) versus \V,, | 
Vis the specific volume at the pressure p at 


nates are (V/V) 
temperature 7, Vy is the specific volume at 
and 0° C, and F is the gas constant in approy 
units. The 
calculated from eq 1. 
to errors in pressure is variable, with the «i 
becoming infinite at zero density. The effi 
117° C of a 0.2 percent error in the press 
illustrated by the distance // between tly 
the chart and the adjacent dotted curv: 
The 
oxygen are represented by the various px 
figure 1. These data are given in Int 
Critical Tables, [3]. They comprise mea 
from three sources namely: (1) Holborn at 


continuous curves represent 


The sensitivity of th 


available measured values of p!/ 


rl 
i 


at the Phy sikalische Technische Reichs 


(2) Kamerlingh Onnes and his ass 
Leiden [3, 5], Amagat [6]. 
(1) The first 


ments at pressures up to 100 atm at 


group consists of nine 1 
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100° CC 
these measurements 


0°, 50°, and 


Cs 


The points 


represent in 
lhe empirical equation represents these 
nts with average deviations from the 
pressures of 0.011, 0.043, and 0.017 
spectively. 

second group are less precise and prob- 
ecurate, Except for those observations 
17° C, they are represented in figure 1 
irious kinds of circles. Table 1 


ire range covered by these measurements 


shows 


the accuracy with which the empirical 
The 


spectively the temperature in degrees C, 


represents them. four columns 


imber of measurements at that temperature, 





the pressure range covered in atmospheres, and 
the average difference between the calculated and 
observed pressures in percent. 

Except for the data at —102.46°C, the average 
deviations at low temperatures do not exceed those 
in the temperature range 0° to 20° C, where the 
equation is more accurately determined by the 

Reference to 
102.46 deviate 
There is no trend 


measurements at the Reichsanstalt. 
figure 1 shows that the data at 
consistently from the equation. 
in the deviations at other temperatures to indicate 
that this is a failure of the equation. It appears 
rather that the data along this particular isotherm 
are in error. Likewise the four observed values at 
113.97° C and the respective approximate densi- 





Comparison of observed 


rcles, Values 


of the chart indicates 


ed by Holbert 
ht H between 


Various ¢ 


ind toy 
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the eff 


al ies 


observed by Onne 


t 
« 


} 


calculated fram equation I 


observed by Amagat [¢ 


houndery 


continuous curves ca 


of field covered by t 





Comparison of values from empirical equation 


with measurements at Leiden 


Average 
deviation 


empera Number Pressure 
ture Theasure rane 


* One measured vy excluded. If this value is included the 


iverage deviation becomes ( 


Four measured v excluded, see text 


ties of 167, 202, 274, and 334, which, according 
to Leiden Communication 169a, form a complete 
series of observations (series XXXII), appear to 
be even more in error. The relative size of the 
deviations are such that the error could either be 
in the temperature or in the calibration of the 
small end of the piezometer. 

(3) The third group of observed data are repre- 
sented in figure | by the letter A. A comparison 
of these data with values calculated from eq 1 are 
given in table 2. No estimate of their accuracy 
can be deduced from their precision, since the 
values given by Amagat were read from smooth 
curves through the plotted experimental data 
However at 100 atm, the values calculated from 
eq | agree with the first groups of data and there- 
At this 


pressure, the deviation from the Amagat data is 


fore must be very near the correct value 


0.3 to 0.9 percent. If the Amagat data do not 


contain errors that increase progressively with the 


pressure, it seems likely that the pressures calcu- 


lated from the equation are within | or 2 percent 
values up to 1,000 atm in the 
temperature range 0° to 200° C 

There appear to have been no new data on the 


of the correct 


critical constants for oxygen since the review by 
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TABLE 2 Comparison of values from eq 1 witl 


ments of Amagat 


The body of the table contains the percentage by wl 


pressure exceeds the calculated value 


remperature, °¢ 


Pickering [7]. The critical density as reporte 
The critical temper 


Dorsn 


about 301 Amagat units. 
ture and pressure reported by Onnes, 
and Holst [8], whose values Pickering chooses 
154.27° K and 49.71 atm 


pressure calculated from eq 1 at this density 


respectively 


temperature is 48.3 atm., whereas at the s: 
density the calculated temperature correspol! 
to49.71latm. is 154.89° K, and the temperatur 
which dp/dV'=0 is 155.27° K. This indicates 
higher critical temperature than that reporte: 
Onnes, although too much dependence shoul 
be placed on slopes calculated from eq 
critical region since, with the number of 
used, the equation is not applicable at densit 
much above those plotted in figure | 

The measurement by Rossini and Frandset 
of the change in internal energy of oxygen bet 
10 and 0 atm at 28° C gives an excellent che: 
eq 1. Those authors report 40 (6.51) 
j/mole. Through the use of eq 1 and int 
of the (OL/OV) p= T(op/oT 
value 258.8 int j/mole is obtained, which ts 


relation 


within the accuracy claimed for the experi 
data. 

K. J. Workman [10] has measured in : 
eter the ratio of the specific heat of o 
various pressures up to 130 kg/cm * to tl 
heat at | atm. <A comparative method 


in which the same stream of gas was pus 
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the calorimeter. Measurements were 

wo temperatures, 26° and 60° C, respec- 
‘hrough differentiation of eq 1 and as- 

that ('p for 1 atm is 7.05 ecal/° C mole, 

of this ratio was calculated for the 

0 kg/em. The third and higher virials 
ected; at the higher pressures it might be 

to include calculations based on the third 

Che values for the ratio calculated for 26° 

i0° C are, respectively, 1.102 and 1.076, 
ereas Workman gives 1.088 and 1.063, respec- 
ev. The cause for this discrepancy of over 1 
cent is suggested by Workman’s statement |10, 


52] as follows: “For example, we can be 
sonably certain that in the present apparatus 
net heat transfer g arises from an excessive 
+ conduction along the line leading from the 
pressure bath to the heat interchanger.”’ 
t transfer at this place would always be in the 
ection to yield a low value for the ratio of 
fic heats 
Equation 1 leads to a value of 94.9 107° for the 
nd virial coefficient for oxygen at 0° C. This 
very good agreement with C. S. Cragoe’s [11] 

timate (95.1 +0.9)* 107°. The value of pV /RT 
0° C and 1 atm derived from eq | is therefore 
4.9 10°°, or 0.999051. Cragoe has shown 
s value to be consistent with measurements at 

esures of | atm or less. The use of R7, 

» 4140 liter atm per mole as given by him leads 
weight of 1.42904 g for the normal liter. This 
i¢ was chosen by Pickering [12] in his review of 
experimental data on the weight of the normal 


Equation | has been compared with the Beattie 
lveman equation [13]. The agreement is very 
. For unit Amagat density in the tempera- 
range 140° to 200° C, the deviation for 
RT at the two extreme temperatures is 8 
ts in 100,000, this being the maximum devia- 
\t 0° C, where eq 1 is in agreement with 
goe’s estimated value [11], the Beattie Bridge- 
n equation gives a value of p\ RT, which is 5 
ts in 100,000 smaller. At higher densities, eq 
n slightly better agreement with the experi- 
ntal data 
It n be of interest to note that the value of 
RT caleulated from eq | for saturated oxygen 
or at the normal boiling point (—183.0° C) is 
which is probably as good a value as any 


, al ttlable 
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III. Virial Coefficients 


Instead of a complete table of pressures over 
the whole range — 200° to +-200° C, for which eq 
1 might be used, values of the second virial 
coefficient B are given in table 3 for every 10 deg 
C, and values of the higher virial coefficients C 
and J) are given graphically for the same tempera- 
ture range in figures 2 and 3, respectively. These 
values are for use with the density in Amagat 
units. If the density in grams per liter is desired, 








200 
TEMPERATURE °C 


s of the third virial coefficuent C plotted as 


function of temperature 


ed with densitic n Amagat units 











TEMPERATURE °C 


ial coe ficient D plotted 


inction of te mm pe ratureé 


ised with densities in Amagat unit 











the density in Amagat units should be multiplied 
by the density at 0° C and 1 atm, 1.42904 g/liter. 


TABLE 3. Values of the second virial coefficient B, multi- 
plied by 1,000, for use in eg 1 


Quadratic interpolation should be used 


remperature, ° ¢ 


Pemperature 


0 10 a» i 40 
Cc 
100 O44 s. 4514 SUH) 4. 4123 5. 0140 
0 0. 488 1. 0743 1, 2125 1. 3647 1. 5333 
1 OSs O. SM45 0.7200 0.6340 0. 461 
10) 140 uO ma + 0027 + (M47 
aw + 2488 
”) ta 7 ™) “) 
loo 5. 7225 6. 048 744 8. 8112 10 AA 
, 1. 7208 1. 9200 2. lt42 2. 4278 2.7247 
) 0. 46%) 0. 3902 0. 3210 0. 250 0 74 
100 0842 1213 1561 INNS 2197 


The quantity pV/RT is dimensionless, and for 
the left-hand side of the equation, p, V, and T 
may be in any units as long as the proper value 
of Ris used. If p is in atmospheres, V in Amagat 
units, and 7 ('+- 273.16, the value of R is 
1/(0.999051) (273.16), or 0.00366434. 

At low pressures the third and fourth terms are 
negligibly small, and when the temperature and a 
pressure not much over 1 atm are given, p/RT 
may be substituted for 1/V in the second term, and 


pV/RT=1—Bp/RT, (2) 


becomes a convenient approximation to eq 1 at 
such low pressures. 


IV. Table and Chart in Engineering Units 


In the compressed gas industry, it is customary 
to express pressures in pounds per square inch, 
temperatures in degrees Fahrenheit, and densities 
in units similar to Amagat units, i. e., the ratio of 
the density at the given condition to the density 
at 70° F and 1 atm (14.696 Ib/in.2). When eq 1 is 
reduced to these units, the constants represented 


by letters are 


B, 1.74509 B,= 0.0009 167 
C"; = 1004.37 & 10° C,= 1.2280 10~ 
dD, 1513.68 * 10° D,=4.1160 10 
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b=150.747/T T= °F +-45:).688 
c= 18.163 < 10°/ T° R—0.0277' 
d= 2636 10/7" 


Table 4 contains pressures that were « 


with the use of these constants in eq 1. mp 
atures are given at the ends of the rows; densi; 
are given at the heads of the columns, and prs 
sures in the body of the table. The uni are 
described in the preceding paragraph excep 
the third section of the table where a unit of 
lb/in.* is used to avoid carrying pressures to 
unnecessary number of places. An attempt | 
been made to give a number of decimal pla 


comparable with the accuracy of the experime: 
data, although some of the values may be ¢ 
to an extra place. 

If the number of pounds mass of gas per cu 
foot is desired, the densities in the units yg 
should be multiplied by the density at 70° F 
| atm, namely 0.082787 lb/ft’. 
duced from three data, namely: (1) th 
0.99933 for pV/RT at 70° F and 1 atm caleulat 
from eq 1, (2) the value 22.4140 liter atm peg 
scale assumed for RT at the ice point, 


Chis factor is 


value 273.16° K assumed for the temperatur 
the ice point. 

The range covered by this table and the ext 
to which the equation is supported by exp 
mental data in this region is illustrated in figu 
by an area bounded on the left by the left-ha 
boundary of the chart, on the right by a ver 
dot-dash straight line, and above and below 
dot-dash curves. 

A chart (fig. 4) has also been prepared that 
produces the data given in table 4 in a form 1! 
permits the rapid solution of practical proble: 
This chart is being printed in a larger size, toget 
with similar charts for nitrogen and hydrog 
and with instructions for their use, as Bu 
Miscellaneous Publication M191. The coo 
nates of the chart are Z pv RT versus press 

On this chart, figure 4, pressures are expres 
in pounds per square inch absolute and temp 
tures in degrees Fahrenheit. Densities ar 
pressed as the ratio of the density at the spe 
condition to that calculated from the equatu 
70° F and 1 atm., i. e., 0.082787 Ib/ft®. Reading 
from the chart indicate that it is drawn | 
sent the equation within 5 Ib/in.? and 
0.0005 for the value of pv RT. The qua 
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SS PABLE 4. Absolute pressures Jor oxygen at various te m peratures and densities 


Densities referred to the density at 70° F and | atm as unity; density at F and 1 atm=0.082787 Ib/ft 


UNIT OF PRESSURE, 1 PSI ABSOLUTE FOR DENSITIES (1 1 
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TABLE 4 Absolute pressures for oxygen at various temperatures and densities Continued 


UNIT OF PRESSURE, 10 PSI ABSOLUTE, FOR DENSITIES (1 \ 
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" aM 215.1 223.9 Zz 242 252 261 21 282 a 
( 224.8 23 244 254 4 274 »% me 3008 
a» 224. ¢ 244. 5 244.6 255 2H 27t 2s7 Jus 10 22 
0 233. & 244.2 254.8 ie V7 Ns 0 112 24 
( “42.8 2 ) 25. 0 276 ss 00 2 2 s 2 
( 252 aK 7 2s7 in) 12 2 ; 2 ‘ 
a ont. ( 273.1 ON5. 4 us i 4 % 9 “ 8] 
0 of 2.7 m5. ¢ tbe) 22 st ( i ist) " 
" 270 202.2 05. 7 lv 4 1s ¢ 7s 4 411 
4 mA wis 1 I w) { | ‘ 4] His 42 { 
i) 206.9 11.3 $25.9 41 sie) 72 {ay 45 422 40 rae 
) ‘ ’ ws ‘ 1 “7 s4 0 41s ht 4 ‘ 
“ 1s “2 wo 2 7 " 4! 4 440 Hi {4s 
" 3.7 ’ 6.0 ry ”) 67 42 444 “ ‘ 
iT 2 ¢ 19.1 iH. 0 tA3 1 10 448 7 437 7 
41.4 s 76.0 4 412 51 1 170 44 512 
1a Lu wi7.9 ws5. 4 PA! 143 " is “4 vt { 0. 
uw) “4 77.3 195.9 415 4 4 475 wt Is ‘4 
14 7.9 ah. 7 405.8 425 +46 $e) {XS TT) 532 
“ “moO 41 7 ie 45 17s iu 22 ' ai * 
: | t 
\ 
> 
— : , , 2 r 
pV RT or Z is independent of the units of p, V, and 2,015 lb/in2 abs corresponds to a densit ; 
and 7’ as long as the appropriate value of R is used. 145.5 and to a value of Z or pV/RT of 0.9415. BP 
The value of Z at 70° F and 1 atm is 0.99933; at we follow the line of constant density for ; 
32° F and 1 atm, 0.999051. For the units given down to 60° F, the pressure read directly fro: 
in the chart, the value of R is 0.0277631. RT is chart is 1,960 Ib/in? abs, or 1,945 Ib/in? gag 
then 0.0277631 (° F+-459.69). The use of the the filling pressure at 60° F. For a more accur 
chart is best illustrated by the following example: estimate of the filling pressure, we note that at 
An oxygen cylinder of 1.528 ft® internal volume same point on the chart the value of Z or p\ 2] 
is designed to be filled at 70° F to a gage pressure 0.9335. We now have at 70° F, pV RI g 
of 2,000 Ib/in.’. 20151,/R(70 + 459.69) =0.9415, and at 60 
a. To what pressure may it be filled at 60° F? PoV2/R(60+ 459.69) =0.9335. Since V,—1\), th 
b. If properly filled, what will the pressure be if leads to p, 72/p2T;=0.9415/0.9335 orp. 2) 
the temperature of the tank should rise to 130° F? (519.69)0.9335 /529.69 (0.9415) = 1960.2. Ibu 
ce. How many cubic feet of oxygen under the or 1945.2 lb/in? gage The direct reading 
standard condition of 1 atm and 70° F, will the the chart is in this case in agreement will 
evlinder deliver, and what weight will be delivered? more accurate calculated value, although ar 8 
d. At 70° F what will the pressure be after in the direct reading as large as 5 Ib/in 
100 {t® have been delivered? expected 
lati (b) If we follow the line of constant densi! 
oe 145.5 up to 130° F, the pressure read directly 
(a) For the purpose of reading on the chart, we the chart is 2,330 Ib/in? abs 
may assume that 2,000 lb/in. gage is equivalent to (c) The number of cubic feet under sta 
2,015 lb/in? abs. The point on the chart at 70° F conditions in the tank =145.5(1.528 2 
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The number of cubic feet delivered is 222.3— Leiden Comm. No. 154a (1918); H. K 
1.528 =220.8 ft®. The weight delivered is 220.8 > C, A. Crommelin, Proe. Roy. Acad. Sei 
0.082787 18.28 Ib 18, 515 (1915); h van Urk « P. G. Nij 
: , . ith Int. Cong. Refrig. 1, 65a and 73a 
(d) After 100 ft.° have been delivered, the Leiden Comm. No. 169a and 169e (1924 
density in the tank is 122.3/1.528=80.0. At this [6] E. H. Amagat, Ann. chim. phys. 29, 68 (189 
density and 70° F, the pressure read directly from [7] 8. F. Pickering, J. Phys. Chem. 28, 97 (192 
the chart is 1,130 |b/in/ abs. [8] H. K. Onnes, C, Dorsman, and G. Holst, P 
Acad. Sci. Amsterdam 17, 950 (1915 
V. References Comm. No. 145b 
[9] Frederick D. Rossini and Mikkel Fra: 
1} R. W. Millar and John D. Sullivan, Bur. Mines Tech Research 9, 733 (1932) RP 503 
Pap. 424 (1928 110] E. J. Workman, Phys. Rev. 332, 1,345 (193 
R. B. Seott, Cyril H. Meyers, Robert D. Rands, Jr., [ll] Carl S. Cragoe, J. Research NBS 26, 4: 
Ferdinand G. Brickwedde, and Norman Bekkedahl, RP1393 
J. Research NBS 35, 39 (1945) RP1661. [12] S. F. Pickering, BS Sci. Pap. 21, 141 (192¢ 
International Critical Tables, 3,8 and 9 (McGraw-Hill [13] J. A. Beattie and O. C. Bridgeman, J. Am. ‘ 
took Co. Ine., New York, N. Y., 1926). 50, 3,136 (1928). 
L. Holborn and J. Otto, Z. Physik. 10, 367 (1922). 
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Reduction of Sphero-Chromatic Aberration in 
Catadioptric Systems. 


By Robert E. Stephens 


Schmidt and Maksutov catadioptriec systems suffer from sphero-chromatice aberrations 


that are of opposite sign. 


By the use of a hybrid system, whose corrector is a deep-meniscus 


lens with an aspheric curve on the external surface, the sphero-chromatic residuals have 


been reduced and, as a consequence, usuable relative apertures have been increased, 


ome time ago a Maksutov camera system of 


relative aperture was designed at the 


onal Bureau of Standards for the purpose of 


rtaining whether or not such a system, em- 
ing as it does only spherical surfaces, could be 
to replace a Schmidt system for certain 
poses. It was found that the residual zonal 
rical aberration was so large that the system 
unsatisfactory. On further analysis, how- 
it was found that the residual sphero-chro- 

aberration of the Maksutov system was 
site In sign to that of a comparable Schmidt 
This 


wested that a hybrid system, partially cor- 


tem and of about the same magnitude. 


d by an aspheric surface and partially by a 
» meniscus lens, would suffer substantially 

‘from sphero-chromatic aberration than either 
Schmidt or the Maksutov type. Such has 
ed been found to be the case. 

\ hybrid Maksutov-Schmidt system consisting 
orrecting plate, plane on one side, placed with 
ertex at the center of the concave mirror, and 
neniscus achromatic lens properly placed, 

d presumably be designed that would be as 

from coma as the best of either of the original 

s. This has not been attempted by the 


hor. Instead, several systems have been 
wgned in which the aspheric surface is applied 
the side of the meniscus lens away from the 
ol 


} 


the chromatic aberration of either a Schmidt 
Maksutov system results from the variation 


orreetion for spherical aberration with index 


k has been sponsored by the Bureau of Ships, Navy Department 


eduction of Sphero-Chromatic Aberration 


of refraction of the correcting plate or meniscus 
lens. In either of these systems the variation of 
spherical correction is | preponderantly due to 
changes of the third order spherical aberration of 
the correcting plate or lens. In a Schmidt system 
the speherical aberration is perfectly corrected for 
one wavelength. At other wavelengths there 
exists simple spherical aberration, undercorrected 
overcorrected for 
Maksutov 


Superimposed upon the 


for longer wavelengths and 


shorter wavelengths. In a system 
the situation is reversed. 
zonal spherical aberration, which is the best that 
can be accomplished by the meniscus lens, is 
overcorrection of spherical aberration at longer 
wavelengths and under correction at shorter wave- 
lengths. Because of the large angles of incidence 
existing at the surfaces of the lens in this type 
system, aberrations of higher order are present to 
a greater extent than in the Schmidt system, but 
the differences over a small range of wavelengths 
are still predominantly of the third order. 

For the first trial it was intended to design a 
system with a relative aperture of (0.7, which 
would be achromatic at //1.0. However, because 
of the effect of the corrector upon the equivalent 
focal length, the achromatic zone turned out to 
be at f 0.89. 


ness were assumed for the corrector. 


A first radius and a suitable thick- 
The second 
radius was then given the value that resulted in 
achromatism at the selected zone. Rays were 
then traced through the system at the most im- 
portant wave length and the condition of the 
spherical aberration determined. Different values 
were assumed for the first radius until a system 
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of the Maksutov type was found in which the 
marginal longitudinal spherical aberration was 
about half that for the uncorrected mirror. 

The first surface was then replaced by an aspher- 
¢ surface * such that the optical path lengths for 
the principal wavelength from all points of an 
incident plane wave to the chosen focus were made 
equal. The chosen focus was that corresponding 
to the ray of the principal wavelength through the 
achromatic zone of the tentative spherical surface. 
The resulting aspheric surface is tangent to the 
tentative surface at the achromatic zone. 

The system was then analyzed for residual 
chromatic aberration by tracing rays at the ex- 
treme waverengths through the system. The 
longitudinal displacements, from the focus of the 
rays of the most important wavelength, are shown 
graphically as a function of height of incidence in 
figure 2. This should be compared with a similar 
graph for a Schmidt system of somewhat smaller 
relative aperture (f/0.83) shown in figure 1. The 
improvement is obviously considerable in spite 
of the larger relative aperture of the hybrid system. 


105 
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Zone Radius in mm 





° 
4S’ in mm 
aberration curves fo a NSchmidt 
amera system 


focal length =83.56 mm 


Although the improvement over previous sys- 
tems is considerable, the balance of the residual 
spherochromatic aberration is not optimum. 


This optimum condition would be attained if the 


s not the purpose of t paper » describe the methods used to con 
aspheric urfaces e pro . sel were imilar to the exact 
wod of F. A. Lue pt. So m. 30, 251 (1940); and, for finite object 


H. S. Friedmar pt. Si m., 37, 480 (1947 
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FicgureE 2 Spherical aberration of first Schmidt-M 


hybrid 


Equivalent focal length= 88.68 mm 


best focus (in the presence of residual aberrat 
for each of the extreme wavelengths coul 
made to coincide with the perfect focus of 
principal wavelength. If this is impossibl 
next best adjustment would be such that the f 
ranges * for the extreme wavelengths over 
Accordingly, further trials were made 

In succeeding designs the second radius of 
corrector was so chosen that modest overcorree' 
of chromatic aberration was obtained at 
selected zone. Then, when the aspheric sur! 
had been substituted for the tentative first surf 
correction of the type illustrated in figur 
sulted. The residual optical path difference 
the focus of the principal wavelength, wer 
than one-tenth wavelength for the extrem 
lengths. The coma of the system, as indi 
by its maximum offense against the sine the 
of 0.007, is greater than would usually be | 
but was quite satisfactory for the purpo 
hand. Figure 4 shows a drawing of this sys! 

The systems designed by the author 
for use in the infra-red region between Su 
10.54 being chosen as the principal way 

4 discussion of the “focal range” is given in Glazebrook 


applied physics, in the section on the opties of the micros¢ 


Conrady), 1V, 221 
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i The principles involved in the reduction of the 
7 aa * ™ residual spherochromatic aberration should, how- 
ever, apply equally well to systems for use in the 
visible or any other spectral range. 


i 





Spherical aberration of Pat hmidt-Maksutor 


hich almost optimum halance of residual 


ymatic aberration has been achieved 


Equivalent focal length= 


=84.66 mm 


WASHINGTON, December 19, 1948. 


eduction of Sphero-Chromatic Aberration 








y.$. Department of Commerce 
Nationa! Bureau of Standards 


Research Paper RP1893 
Volume 40, June 1948 


Part of the Journal of Research of the National Bureau of Standards 





Some Energy Relations in the Systems PbO-B2O3 and 
PbO-SiO2 


By Leo Shartsis and Edwin S. Newman 


Heats of solution of a series of glasses in the svstem PbO-B,O; were measured in 2 N 


nitric acid with a simple vacuum-bottle calorimeter. 


Heats of solution of a series of glasses 


in the system PbO-SiQ», ranging in composition from 6 to 35 percent of SiO», were measured 


> 


in a mixture of HF and 2.5 N nitric acid with an electrically calibrated isothermal-jacket 


calorimeter 


From these data, and heats of fusion available in the literature, calculations 


were made of the heats of reaction of the component oxides in either amorphous or crystalline 


states in forming glasses in these systems. 


I. Introduction 


heat accompanying the glass-forming re- 


ns is of interest to the physical chemist 


dying glass. Its significance lies in the fact 

a» knowledge of the energy changes accom- 
ving glass formation should ultimately lead 
a better understanding of the glassy state. 

heat of glass formation is of practical im- 
tance, because it is one of the factors needed in 
luating the efficiency of glass-melting furnaces. 
lhe direct measurement of the heat accompany- 
sthe reaction of components to form a glass is 
vy difficult, if not impossible, because the re- 
on takes place at elevated temperatures and 
ften very sluggish. For this reason, indirect 
surements were made by use of a heat-of- 
ition method. The reactants were dissolved 
i suitable solvent and the heat of solution 
rmined. The finished glass was dissolved in 
sume solvent in such an amount that its solu- 
was identical with that resulting from the 
tion of the reactants The difference between 
heat of solution of the reactants and that of 
finished glass was the desired heat of reaction. 
the calorimetric procedure was conducted near 
m temperature, the heats of reaction so de- 
mined were at room temperature (approxi- 

25°C). The heats of reaction at elevated 
peratures may be calculated from such data 
bined with the specific heats of the reactants 
product over the desired temperature range. 


Heats of Reaction of Some Oxides 


II. The System PbO-B.O, 


1. Preparation of Glasses 


Two series of glasses were used in this study. 
The first was made specifically for heat-of-solution 
measurements; the second was made for another 
project, and portions remaining from that work 
were used to make additional heat-of-solution 
measurements. The two series were prepared in 
essentially the same manner, but details varied 
somewhat. Only the preparation of the first 
series will be described. 

The raw materials used in melting the glasses 
were boric acid and PbO! of reagent quality. 
Fifty-gram batches were melted in platinum cru- 
cibles heated in an electric muffle to not less than 
800° C. For melts high in B,O,, temperatures as 
high as 950° C had to be used. During melting 
the glasses were homogenized by frequent stirring 
with a platinum wire mounted at the end of a long 
glass tube. For very fluid melts (those high in 
PbO), the mixing was aided by shaking the cruci- 


When 


it appeared homogeneous the melt was quickly 


ble and moving it with a swirling motion. 
poured onto a clean steel plate. Each glass was 
ground to pass a No. 60 sieve, well-mixed and 
remelted, 

Glasses containing more than 92 percent of 
PbO could not be obtained because of crystalliza- 


ray pattern made by H. F. MeMurdie indicated that t 
te, the 8 or low-temperature form 








tion. Geller and Bunting [1]* obtained clear 
glasses containing as much as 94.5 percent of 
PbO, but they prepared much smaller quantities 
of material and used a more efficient quenching 


system. 


2. Heat Treatment and Chemical Analysis 


Most of the glasses were annealed by heating 
to 20°C or 30°C below their deformation tempera- 
tures [1] for an hour and then cooling very slowly. 
The heat of solution of a strained glass may be 
expected to be higher than that of the same glass 
after annealing. However, the method used was 
not sufficiently sensitive to detect this difference. 
Therefore, the remaining glasses were not an- 
nealed 

The glasses were analyzed for PbO. The PbO 
content was determined as the sulphate by several 
evaporations with HF and H,SO, and heating at 
500° to 550° C to constant weight. 


3. Calorimetric Materials, Apparatus, and Pro- 
cedure 


A quantity of 2.00 N nitric acid was prepared 
from chemically pure acid and distilled water. 
About 500 ¢ of analytical reagent grade ZnO 
was heated for 1 hour at 900° to 950° ¢ 
and ground to pass a No. 100 sieve and stored in a 
This material was 
used to determine the effective heat capacity of 


’, cooled, 
well-stoppered reagent jar. 


the vacuum-bottle calorimeter. 

Most of the 
were made with a simple vacuum-bottle calorim- 
eter [2] 
wide-mouthed Thermos flask in a container of 


heat-of-solution determinations 


As can be seen (fig. 1), it consists of a 
insulating material. Inserted through holes in 
the cork stopper are a two-bladed glass stirrer, a 
Beckmann funnel 
through which the charge is introduced into the 
disk marked with six 
equally spaced broad radial lines was glued to the 


thermometer, and a glass 


calorimeter. A paper 
top of the horizontal member of the friction drive, 
and a 3-neon bulb (not shown) was used to illumi- 
nate it. When the stirrer is operated at a speed 
of 600 rpm, the lines appear stationary. At 


higher speeds the lines appear to move in the same 


Figures in brackets indicate the literature references at the end of this 


paper 
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direction as the stirrer, whereas at speed~ low. 
than 600 rpm they appear to move in the reyers, 
direction. A speed-adjusting screw permits thy 
speed to be maintained approximately cons: ant 
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FIGURE 1 Vacuum-bottle calorimeter I 


As ZnO was used to measure the heat capacit) 


of the calorimeter, its heat of solution in 21) 
N HNOs was determined with an electrically cali- 
brated calorimeter, which will be described lat 
As a check on the vacuum-bottle calorimeter, t! 


heats of solution of crystalline PbO and of 


/ 


PbO-B,O, glasses were also determined with 


apparatus. l 
he 

4. Procedure and Calculations » 

. th 

The acid charge consisted of 425+0.0 * 
2.00 N nitric acid. Its temperature was adjusted e 


to 2 or 3 degrees below room temperature belo! 





introduction into the calorimeter, so that 
initial temperature indicated by the Beekman 
thermometer would be at the lower portion of !ts 
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er the calorimeter was assembled a acid. For subsequent experiments with glasses 











lower sing riod of not less than 10 minutes was al- or oxide mixtures, their heat capacities were sub- 
Vers ved he reading of the Beckmann thermome- stituted for that of the ZnO in the value obtained 
di was en recorded, the clock started, and the for the calibration. 
™ ample vowdered to pass a No. 100 sieve) intro- The solid material introduced into the calo- 
ed a uniformly as possible in not less than 1 rimeter was at room temperature, which was 
ite more than 2 minutes. Temperature higher than the initial temperature of the calo- 
dings were taken every 5 minutes until the rate rimeter. Account was taken of the sensible heat 
remperature change of the calorimeter was thus carried into the apparatus. This heat 
; stant for at least three 5-minute intervals. is given by the formula C, (T,—To), where 
experiments showed that with a properly designed C, is the heat capacity of the solid, 7, is the 
rer and an adequate stirring rate, complete temperature of the room, and 7 > is the initial 
tion Was attained in 20 minutes, so that in temperature of the calorimeter. In the determi- 
ral the total time of the experiment was 40 nation of effective heat capacity, this amount of 
ites, of which the last 20 were used as a rating heat must be added to that developed by the 
" reaction of the zine oxide and acid to obtain the 
bd he caleulation of the final temperature was total heat input into the calorimeter. For the 
by extrapolating the constant rate of tem- heat-of-solution measurements, this amount of 
ture change of the calorimeter during the heat must be deducted from the total heat 
yw period back to zero time. The “final” measured to obtain the quantity solely due to the 
: perature so calculated minus the initial tem- reaction of the sample with the acid. 
ture of the calorimeter was taken as the tem- ; 
iture rise of the calorimeter caused solely by S. Data and Discussion 
energy evolved by the reaction of the sample Table 1 gives the data obtained, and figure 2 
the acid. This corrected temperature rise shows them in graphical form. As the heats of 
“ the calorimeter was converted to an energy solution of the first two mixtures of PbO (crystal) 
nge by multiplying it by the effective heat 
wity of the calorimeter and its contents. The 
rs inherent in such a method of calculation are 
woximately the same in both heat-of-solution 
t rminations and in the effective heat-capacity : 
rminations and thus tend to compensate each a Se 
= : 
* lo determine the effective heat capacity of : 
) calorimeter, a heat-of-solution experiment 
- ZnO as the sample was made. The known ~ 
late t evolved by ZnO [—256.9+0.1(T-—24)] eal/g, = 
etl temperature in deg C—previously determined 
f woth an electrically calibrated calorimeter, when 
» this ded by the observed temperature rise vielded 
leat capacity of the calorimeter. Since the : - 
hermal heat of solution of ZnO varies with PDO, PERCENT 
perature, the value at the final temperature Ficuns 2. Heats of solution in the system PbO-B,O, 
the calibration was used. The effective heat 
v . . Mixtures of PbO (gls) and ByO; ‘gls) (curve 1); mixtures of PbO (c) and 
a neity obtained was thus that of the calo- adie telich dnasein es snietincan al Deis taiehs aide idi t Santee catalina 
ter and the reactants, 1. e., the ZnO and the of PbO (c) and BzOs (¢) (cu-ve 4); PbO-B20s glasses (curve 5). Open circles 
elo! except 0 percent of PbO) represent values obtained with 10-¢ samples 
_i— closed circles represent values obtained with samples of 18 to 25 g. The 
ry defined as 4,1840 abs j is used in this paper. According to difference between the heat of solution of PbO (gis) and PbO (ec), repre 
mann heat evolved is negative sented by F, is the heat of fusion of PbO at room temperature 
of its 
Heats of Reaction of Some Oxides 473 


2arch 








lanie 1. Heats of solution in the system PbO-B,O plus B,O, (glass) tested fell on a stra +) 
between the end members, no other de orm) 


Vacuum-bottle calorimeter | Precision calorimeter tions of oxide mixtures were made, 
— A curve (5, fig. 2) was fitted to the heat. of so 
Sul Percent-| Heat c Number tomate Number , a - ; pa . 
seat | sol. |. Ofde- | onie «| Heat of | of de tion of the first series of glasses. The jjeats , 
POO | tha oy ee a solution of the second series (solid circles |: fig 9 
were somewhat lower than those of the first gop, 
FIRST SERIES because the amounts of sample dissolved jy 4) 
- ee latter were about twice as great as in the forn, 
BOs} o| —ae 1| on thus causing the final concentration to be differey4 
san _ = » The heats of solution of glasses having high per 
De — 6 3| w centages of PbO seemed to be linear with th 
Do 0.9 —43.5 1; weight percent of PbO. The equation of a |iy 
oe op i “s 55.2 ; was fitted by the method of least squares to th 
“ 85.4 65.7 3, 0 data obtained with the vacuum-bottle ealorimer, 
“~ —e- aon - “ on glasses containing more than 85 percent 
De r 79.3 PbO The value at 100 percent of PbO obtai 
Licharetee 100 ( Sr. : 10 87.7 2 . oo 
een idigaa from the equation of this line was —99.5 ¢a| 
-—~ +PbO so. " . . Since the average vacuum-bottle value of the hea 


of solution of PbO (c) was —87.6 ecal/g, the hea 
of fusion of PbO (ce) at room temperature caley 
lated from these data is +-11.9 cal/g. 

The value for the heat of solution of B.O 


SECOND SERIES 


(ila “2 l IS 
De ape v8 2 18 was taken as —53.5 cal/g, the difference bet 
ns pe - : — the heat of solution of glassy BoO, (—116.1 eal y 
we 7.7 12.5 2 2 and the heat of crystallization of B,O, (c) at roon 
~ 70.9 “ - temperature (—62.6 calg [3] ). Straight li 
ae a2. 1 i. 4 : 25 were drawn between the values for the heats 
“4 ~ p : > solution of B,O, (gls) and PbO (gls), curv: 
Do a8. 0 70.2 2 2 B.O (gls) and PbO (ec), curve 2: B.O, (c) an 
oe na m . - os PbO (gis), curve 3; and B,O, (c), and PbO 
Do loo(el W. 5 curve 4. These lines represent the heats of solu 
tion of mixtures of the oxides. These data 
Rpts ew a re gether with the heats of soluiion of the glasses 
2 In 2.00 N HNO enable a calculation to be made of the heats a 
Rivero tangerine ev?“ ea gt alaaae companying the following four reactions at ap 
Extrapolated proximately 25° C: 


rPbO(gls) + (1 — 2) B,0,(gls) =zPbO.(1 —2x)B,0,;(gls) + AH, 
rPbO(e)- (1 r)B.0,(gls) =rPbO. | 1—z)B,0 (gis) ' AH, y 
rPbO(gls) + (1 — 2) B,0,(c) = 2PbO.(1—2z)B,0, gis) +-A//,; 


zPbO(c) + (1—2)B,0,;(c) =2PbO.(1—2z) B,O,(gls) + AM, 
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“ais a 


searca 


chich - and (1—r) are the grams of the appro- 
ite 0) de necessary to make 1 g of glass and 
H,, et the heats of glass formation, are 
viables dependent upon x. These calculations 

mate by subtracting the heats of solution 


‘the glasses from those of the appropriate oxide 
xtures on the assumption that the heat of solu- 
of the mixtures is a linear function of the heats 
‘lution of the component oxides. The results 
shown graphically in figure 3. The two curves 


re PbO (gis) is one of the reactants are hypo- 
tical, because glassy PbO at room temperature 
syet unknown. Table 2 gives values of heats 
action at such intervals that linear interpola- 


mav be used. 


» Heats of glass for mation at room temperatures 


wimately 25 ( from oxides in the system 
-KO 
Reactants 
P PbO (el PHOK PbO (gels PHOK 
BeO (gels B,Os(gls BpOg(« BeOsl« 
age 
” Heat of reaction, cal/g of product 
All Alhk; SH Al, 
0 0 +62. 6 (2.6 
if 6.0 +52.9 M5 
13.1 11.9 43.2 44.4 
1v 17.9 1. 5 
As] 26. 2 23.8 23.9 26 
2.8 29.9 14.1 17.1 
4 16.9 +4.4 0 
1.0 41.58 ; 1.1 
2.7 is. 9 15.1 10 
SS MS 24.4 19.0 
f2 rv’ 10 25.0 
t me th. 5 70 
t ri “fe ik. 7 16 
hie 41.3 t 
61.2 2.9 2.4 4.1 
4.9 17.9 41.2 r4 
™ is 5 tA. 8 ri 
“.¥ i! 7 Is. 4 
25.4 14.6 19.1 a4 
12.7 1.4 of 0.7 
( 11.9 0 11.9 
40 abs 


convention, heat evolved is negative 
rd, J. Am, Chem, Soe. 68, 3147 (1941 


of Reaction of Some Oxides 




















-60 T T T T 
-60} : 
| 
“40 4 
= 
<a 
« 
rT) 
« 
Y -20 7 + 
o , 
aw * 
« J . 
° ‘ ‘ 
_ ‘ ‘ 
Ps} Oo, i 
z PbO (ga) +B O fc) ,’ | 
4 23 4; ‘ 
: ON 
ty 
! .eoL ha Peo tiv, 0. @) J 
a 
~ Ps 
< f 
w sf 
=z 4’ 
+40 F - 7 
i! 
$ 
/ 
f 
/ 
+60 }* 4 
1 rl 1 1 
° 20 40 60 80 100 
PbO, PERCENT 
Figure 3. Heats of reaction to form glasses in the system 


PbO-BoOs. 


Ill. The System PbO-SiO, 
1. Preparation of Glasses 


The glasses investigated were prepared by melt- 
ing, in platinum crucibles, the requisite amounts of 
PbO and SiO,. The weights of the melts ranged 
between 300 and 500 g. They were stirred with 
a platinum rod and when they seemed homo- 
geneous were poured into iron molds. Glasses 
containing more than 35 percent of SiO, could 
not be conveniently melted in the apparatus 
available. Melts containing less than 13 percent 
of SiO, crystallized when poured into molds. 
In order to prepare glasses containing less than 13 
percent of silica, small portions of the crystallized 
compositions were melted in platinum crucibles 
and poured into ice water. In this manner glasses 
containing as little as 6 percent of SiO, were 
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prepared. Microscopic examination of the lowest 


silica-containing glass thus quenched showed it to 
contain less than | percent of crystalline material. 


2. Heat Treatment and Analysis 


The glasses were given no further heat treat- 
ment. They were analyzed for PbO in the same 
manner as the PbO-B,O, glasses. Silica was ob- 
tained by difference. 


3. Calorimetric Materials, Apparatus, and Pro- 
cedure 


A quantity of 2.50 N nitric acid was prepared 
from chemically pure acid and distilled water. 
Reagent quality hydrofluoric acid (about 48 per- 
cent of HF) was used. The calorimetric charge 
consisted of 50 ml of HF and enough nitric acid 
to make 648.5+0.05 ¢ of mixed acids. 

The calorimeter (fig. 4) used in this work has 
been previously described [4]. The whole assem- 
bly was placed in a water bath whose temperature 
was kept constant within +0.005° C by auto- 
matic means The temperature of the calorim- 
eter and water bath were measured with 
platinum-encased platinum resistance thermom- 
eters whose resistances were determined with 
Mueller temperature bridges. 

The Regnault-Pfaundler formula [5] was used 
in calculating the temperature rise obtained and 
to guard against arithmetical errors, the tempera- 
ture rise was also calculated by the second Geo- 
physical Laboratory formula [5], by using the 
value of the thermal-leakage modulus obtained by 
the former method. 

The effective heat capacity of the calorimeter 
and contents was measured by introducing a 
measured amount of electrical energy and meas- 
uring the resultant temperature rise. 


4. Data and Discussion 


Table 3 gives the data obtained, and figure 5 
shows them in graphical form. The equation of 
a straight line was fitted by the method of least 
squares to the values of heats of solution of com- 
positions containing less than 17 percent of silica. 
The value at 0 percent of silica (100 percent of 
PbO) obtained from the equation of this line was 

99.2 calyg. As the value for the heat of solu- 
tion of PbO (c) was 89.4 cal/g, the heat of 
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Heat of solution in the system PbO-Sif do. 

PbO (uls (curve 1); PbO-SiO, glasses (curve 2 
rhe difference heat of 
represented by F, is the 


ind SiOz» 
{PbO « ind SiOz between the 
PbO (gis ind PbO (« 


on of PbO at room temperature 


curve 3), 


pted. The heats of solution of silica gel and 
ous silica were determined in a high-HF con- 


a 
The 


corrected 


difference between these two 
for the heat of dilution 


sed by the water contained in the gel the is 


tration. 
es when 


ey difference between these two forms of silica 
therefore independent of the solution in which 
as measured. This difference was then added 
wbraically to the heat of solution of silica gel 
i lower concentration of HF, such heat being 
lated to the ignited weight of the gel and 
ected for the heat of dilution caused by the 
The necessary heats of dilution 
The heat 
wlution of quartz was calculated by subtract- 

36.8 calories (the value obtained by Mulert 


in the gel. 


measured in separate experiments. 


rPbO(gls) + (1 —x)SiO,(gls) 


rPbO (c)+ (1—2) SiO, (gls) 
rPbO (gis) + (1 —2r) SiO, (ce) 


rPbO (ec) + (1—7) S10, (ce) 


of Reaction of Some Oxides 


rPbO.(1 


[6] for the difference in heat content of amorphous 
silica and quartz at room temperature) from the 


heat of solution of vitreous silica. 


TABLE 3 Heats of solution in the system PbO-SiO, 


Approxi 
mate 
weight of 
sample 


Number 
of deter 
minations 


Heat of 
solution 


Percentage 


Substance of SiOz 


4) ml HF +sufficient 2.5 (N) HNOs to make 648.5 g 


cal/g 
PbO (« 
PbO (gis 
(ilasses 
Do 
Do 
Do 
Do 
Do 
Do 
Do 
Do 
Do 
Do 
Do 
Do 
Do 
Silica gel 
H,O 
SiO» (gls 


SiO, (quartz 


250 ml H F+sufficient 2.5 (NV) HNO, to make 648.5 g 


Silica gel 
SiO», (gis 
H,O 


4.1840 abs j 

? By convention, heat evolved is negative 
Extrapolated 

‘ Calculation made on the basis of ignited weight 

> Calculated 

*0. Mulert, Z 


1 calory 


inorg. Chem. 75, 198 (1912 


With these data and the previously outlined 
of the 
accompanying the following reactions at approxi- 


methods, calculations were made heats 


mately 25° C: 
rPbO.(1—2)SiO,(gls) + AH, 


rPbO.(1—2)Si0,(gls) 


AH, 
rPbO.(1 AH, 


+ AM, 


r)SiO,(gls) 


r)SiO,(gls) 





Ficure 6 


TABLE 4 








Heats of reaction to form glasses in the system 
PbO-SiO, 


Heats of glass formation at room temperatures 
approximately 25° ¢ from oxides in the system PbO-Sit )» 


Reactants 


PbO (gls)-+ PbhOl(ce)4 PbO(els)+ PHOkK 
Product SiO(gl SiOe(gls SiOd(« SiO 
glass 
percentage 
of SiO HEAT OF REACTION, CAL/G OF 
PRODUCT 2 


Mit, 


1 calory = 4.1840 abs j 
According to convent 
O. Mulert Z. anor Chen 


In the 
of data, straight lines were drawn betwee; 


The results are shown in figure 6. 


100 percent of silica. Table 4 gives values 
of reaction at small enough intervals th; 
interpolation may be used. 

By combining the data from both the 
+2.4 keal mole 
tained for the heat of fusion of PbO at room ; 
perature. 


investigated, a value of 


This may be compared with a repor 
value of + 2.75 keal/mole at the melting p: 

It is to be expected that the heat of fusion at ; 
temperature should be lower than at | 

the melting temperature. 


IV. Summary 


The heats of solution of some compositions i) 
systems PbO-B,O, and PbO-SiO, have be 
sured. From these data and from values of 
of fusion available in the literature, heats of » 
tion of appropriate oxides to form glasses in ¢} 
systems have been calculated. 
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Divided Flow, Low-Temperature Humidity 
Test Apparatus 


By Arnold Wexler 


An apparatus is described, based on the principle of divided flow, for producing air of 
known relative humidity at temperatures below 0° C with an average error of 3 percent, 
It was designed primarily for testing and calibrating at low temperature the electric-hygrom- 
eter elements used in radiosondes, although any device with a diameter less than 1\ in. may 
be inserted into the test chamber of the apparatus. It is capable of producing rapid and 


discrete changes in relative humidity at constant temperature, a desirable feature for studv- 


ing lag characteristics. 


1,500 ft/min 


I. Introduction 


\n apparatus has been developed for producing 
f known relative humidity at temperatures 

0° to —40° C. The primary purpose for 
ling this apparatus was to provide a means 
nvestigating the behavior of radiosonde hy- 
ters at low temperatures, especially electric- 
ometer elements of the Dunmore type. In 

use, on radiosondes, the elements are sub- 
d to rapidly changing humidity conditions 
to air velocities of about 700 ft/min. It was 
lesirable in studying and calibrating these 
ents to be able to do so under simulated 
tions of use. Therefore, in the design of 
apparatus, provision was made to permit 
change from one discrete value of relative 
dity to another, and a device was incorpor- 
for controlling the velocity of the air. 


II. Theory 


method utilized for producing an atmos- 
of known relative humidity was to divide 
ately a current of dry air into two streams, 
{ which was maintained dry and the other 
ated with respect to ice, and then to recom- 


the two. It was employed by Gluckauf ? in 


opment was financially supported by the Bureau of Ships, 
{the Navy 
wuf, Proc. Phys. Soc. 59, 344 (1947 


midity Test Apparatus 


It is provided with a control for obtaining any air velocity up to 


an investigation on absorption hygrometers and 
by Walker and Ernst, Jr.6 in a laboratory setup 
for preparing constant mixtures of air and water 
vapor. Figure | is a simplified schematic diagram 
of the basic components of the apparatus illustrat- 
ing the principle of operation. By means of the 
proportioning valve, VV, a flow of moisture-free 
air is divided into two parts in a known ratio. 
One part is passed through the saturator, S, over 
a series of trays containing ice until it is completely 
saturated. It is then mixed in the mixing cham- 
ber, Cy, with the other part that has been main- 


tained dry, and allowed to exhaust through the 


ATMOSPHERE 


{ 








TEST 
CHAMBER 


Cy 


PROPORTIONING VALVE 
VY 

















= | 
ORY MIXING 
AIR CHAMBER Cm 





CONSTANT TEMPERATURE BATH 


FiGuRE 1 Simplified schematte diagram illustrating the 


princi ple of operation of the humidity apparatus 


C. Walker and E. J. Ernst, Jr., Ind. & Eng. Chem. Anal. Ed. 2, 134 
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test chamber, C;, into the atmosphere. The 
saturator, mixing chamber, and test chamber are 
kept immersed in a constant temperature bath. 
In order to subject a hygrometer or other device 
to air of known and constant relative humidity, 
it is inserted into the test chamber, C7. 

The relative humidity in the test chamber is a 
function of (a) the fraction of air that passes 
through the saturator, (b) the total pressure in 
the saturator, (c) the total pressure in the test 
chamber, (d) the saturation pressure, and (e) 
the partial pressure of the water vapor in the test 
Under conditions, which the 


chamber. ideal 


apparatus approaches with sufficient closeness, 
all variables except one drop out, and the relative 
humidity becomes equal to the fraction of air 
fhat passes through the saturator. If Y is the 
fraction of air that passes through the saturator, 
HW’, is the ratio by weight of the water vapor to 
the dry air of the mixture leaving the saturator, 
and W, is the ratio by weight of the water vapor 
to the dry air of the mixture passing through the 
test chamber, then 


W-= XWs. (I 


From Dalton’s law of partial pressures, and 
assuming the perfect gas laws apply, it follows 
that for an air-water vapor mixture the ratio by 
weight W of the water vapor to the weight of dry 
air is 

M, p 


WW M, P p’ (2) 


~ 


where Mf, and M, are the molecular weights 
respectively of air and water vapor, p the partial 
pressure of the water vapor in the mixture and P 
the total pressure of the mixture. Substituting 


the value of W in eq 2 into eq 1 yields 


Y (P. Ps)Pe. = 
: (Py Pc) Ps 

Relative humidity is defined as the ratio of the 
actual vapor pressure to the saturation pressure 
of water. The relative humidity, //, in the test 


chamber, C,, is therefore 


Haut. (4) 
Ps 


480 


Combining eq 3 and 4 leads to 


. 





X 
Ps Pps fo 
( P. P- ) Po— pr 


If the pressure drops from the saturator to 


H 


a 


test chamber and from there to the atmosphe 
are made small, then Pg= P-=P,, the atmosph 
pressure. At 0° C, ps/Poe= ps/P4=0.006, at 
it equals 0.0025, and at —20° C 
Also P-/Pe— po=P./Pa—Pe, 

is less than ps. Therefore, for the applicabl col 


() 
it equals 0.0 
where pe equals @ 


ditions eq 5 reduces without significant error tol 


H=X. 


III. Description of Apparatus 


The essential functional units of the appara 
are (a) the drying system, (b) the proportionii 
system, (c) the humidifying system, (d) the mixif 
chamber. (e) the test chamber, (f) the cooling S| 
tem, and (g) the thermo-regulating system 


Figure 2 


FS 


temperature control. schematica 


shows these functional units and their compondll 
parts. | 
Air is supplied to the drying system at pressuf 
from 70 to 90 psi. Water from the air at he 
supply pressure is removed by freezing in the 7 
ice-varsol bath, C7, in which the temperatur 3 
well below —70° C. The vapor pressure of ’ 
water in the air emerging from the drying syst n 


upon expansion to | atmosphere, is probably 
than 0.0003 mm Hg. Particles of dirt, ice, oF 
are caught by the filter, A8. 

The dry air is brought to room temperature® 
means of the electric heater, A/0, is reduced a 
controlled in pressure by regulator, A/2, and ([ 
passes through the proportioning system. * 
pressure regulator also serves to control the sp 
of air in the test chamber, A245, permitting 
justments up to 1,500 ft/min. 

The proportioning valve, A/4, figure 2. wii 


i 


divides the air in a definite ratio consis 

orifices of equal cross-sectional area so arralt 
that by a turn of the knob of the valve. th 
coming air can be divided to produce any of * 


ratios, 0, 1/6, 1/3, 1/2, 2/3, 5/6, and TI 

ratios are the fractions of air entering \ 

that emerge through one of the exit cl . 
. 


Journal of F ese Hurnid 











69 Aseijog 





O19 Aojey 


snposvddp fipprany yu) fo 














9 
© os40 esnjosedwe uoysuo So ome 
WO 10K ——--——— 2 SS ae en ae ---- JOKOW =—4009 
= JOS40A 
e | 
. | 
= oz¥ l@v [ 3 e2y equious 1 
—— 49|00904g =: 4050884 js > buxiw | | 19 weg 02) Aug 
3 3 a orl Dee 
$ ln y | | [ © | @2) fag puo Leotee | 
$ 80 je: ; | | "h 
° povenBosousees is e si l l | a 
a° =|— 
= gi— PID 81109 
1s ; = | | it | 
1044145 jOBny sued > > - } - | | JOS4DA | - ° 
uw eo $ | l , al 
LJ -} : 
99 + s | | Cceeaataeed 
4944S 104;0004g 9 = ae Ee - a | 
“ = €id "og pies | 
rs) 
W3LSAS z| 12 = 
ONILVINOSYOWNSHL 3] |2 Loa “a 
| e| |e | = 
La ee ae ea cm me ae ee ae ee ae ne a ce a ce me > an ee ee ae ee eed 
19 weg 02; fag 
sai miei ee ee 
° 5 
Buiuoipsodorg Ziv iy = A a > | o ae 
r 7 « > wy dou 
soy ojnBey — /2#01N Bes Oly hy, ain ws tv ugg © 9 Sy © o2nisiow | 2 S 
eunsseig owed 1 194004 | 9 0 me rl O 4109 seddoy >S 5 
6iv Sa fl papane lo © am. ¢ i 
SOALOA UO!JOZIVINDZ sesnssesg {#5- iY T nl | T = 
| 8 + = 
v Vv =e a2 
| i U 0 
Liv @A1DA @)peeN cry basa ac | | bs a < ° t “a @ 2 
p "in F ee 
W31SAS SONINOILYOdOdd oiy srowounn onal } ° os =$ | 
Aunduew L due s 9 
viy e609 | 6y e605 Sy e609 , . €y #605 - ° 
aunssaug ~) eunssasd ( einssard ( ¥/ ) e4NSseig 
W3LSAS ONIANYG JYNLVEYSdW31 MOT 
ee eo Slane — age es, &. BABA SB BS B ae 


All 


sjipd yu ynod Uo pup spun pouortouns y) fo bua pip 


W3iSAS SONIAJICINNH 















VLPDULIYIS Z% aNaol 


W3LSAS ONILVINOWID TOSYVA 






























- + @ 2 hrs 2 > « = tm SS fe 


























Ayddns ay 





Humidity Test Apparatus 


eG 














the valve. If they are considered to apply only to 
that part of the air that subsequently becomes 
saturated, then they may be used as the values of 
NV in eq Sand 6. To assure that the air will divide 
in accordance with one of these riven fixed ratios, 
the downstream pressures must be equalized. 
To achieve this, the differential mercury manome- 
ter, A/6, and the two variable resistances (pres- 
sure equalization valves), A/S and <A/9, are 
installed 

The construction of, and the flow of air through, 
the proportioning valve are shown in figures 3 and 
4. Essentially the valve, A/4, consists of a brass 
base, B, a brass top, 7, 1 steel orifice plate, P, and 
a Bakelite knob, A. The 


valve through a central hole in the base, B, 


incoming air enters the 
con- 


tinues through two 3/16-in. diameter holes in the 


Fiaure 3 ( ect é of the proportioning ure 
ho ng the flo of ai 

Pt rt i rif " K, Bake knob; O 
S \ k nut 





T 


FiGuRE 4, 


R, Brass base; P, ste« rifice f S, steel shaft; T, brass top; V 


422 


Esse niial parts of the proportioning valve 


teel packing ! 
Kt 


midsection of the orifice plate, P, strikes | 


7, reverses its flow and again passes throuh the 
steel orifice plate, P, this time through six |-mm 
diameter orifices, and exits in two streams from 
the base, B. The base, B, has two cireylye 
grooves on a 1's in. diameter, 3/16-in. wide, \; jy § 


ws . 
Each 00 


The orifice plate, P 


deep and separated \ in. at the ends 
is connected to an outlet, O. 
rests on the base, B, in such a fashion that the 
orifices communicate with either or both grooy 

depending upon the angular position of the plat, 


P, with respect to the base, B. The surfaces offi 
contact between the plate, P, and base, My 
ground and lapped and coated with a thin film , 
stopcock grease to provide a leaktight but rotatl 


By means of a steel shaft, S, t] 
connects the orifice plate, P, to the Bakelit: cnol 
K, through the valve top, 7, and a steel packiy 
nut, .V, that seals the shaft, S, in the top, 7, the 
orifice plate, P, 


able joint 


a - 


can be rotated to any one of thell 
seven possible settings Upon turning the knob& 
K, two spring actuated snap pins, Q, in the knob§ 
K, click and seat in a series of holes on the out 
surface of the valve top, T, corresponding to thea 
seven desired positions of the orifice plate, ?, withll 
respect to the valve base, B 

The two air streams, upon leaving the propor- 
tioning system, flow through the humidifying sys 
tem in parallel channels, thermally in contact with : 
one another to allow heat interchange and tem- 
perature equilibrium between the two streams to 
Saturation is ac- 


be achieved and maintained 





ee ee ae 


it; A, Bakelite knob; Q, spring-actuated snap pin; C, steel pin 
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omp shed in four stages. The air to be humidi- 
fed is first brought to ambient temperature in a 
4, A120, figure 2. It then passes through 
iturators, A21, A22, and A23. The pre- 
A20, of staggered shelves, one 
Each 
wries iS pneumatically isolated from, but in ther- 
ntact with, the other. 


reco 
thre 
ree 
oole 


~ries for each of the two currents of air. 


consists 


al The saturators are 
similarly built, except that trays replace one series 
f shelves. The trays can be filled with water 
frozen. Thus, in 


one channel air passes over trays of ice and in 


nd the water the saturator, 
parallel adjoining channel, air passes over 
From the final A23, the 


jurated air and the dry air are centrifugally 


! 
selves. 


saturator, 
xed in the mixing chamber, A2/, and dis- 
arged into the atmosphere through the test 
amber, «125. 
The humidifying system, the mixing chamber, 
and the test chamber, A25, are immersed 
a varsol bath, C2, whose temperature is closely 
rulated. This is done by circulating the varsol 
the bath through cooling coils, C14, immersed 


dry ice. The bimetal thermoregulator, C8, 
ough the relay, C70, controls the heater, C11, 


the bath 


keep termperature at any desired 
w. The bath, C2, is actively agitated. As 
aid in averaging and stabilizing the variations 


the bath the intermediate and 
alsaturators, 422 and A23, the mixing chamber, 


and the test chamber, A295, 


temperature, 


are encased in 
copper shield box, C73, to provide a unit of 
ch heat content that assumes the mean tempera - 
of the bath. 
so] away from direct contact with the enclosed 


The shield box, C13, keeps 


iponents. 

The thermoregulating system is capable of 
trolling the temperature of the humidifying 
to +0.05 is? © -().] 
mm to —30° C + 0.2° down to —40° C 
runs up to 10-hr duration. For short periods 
time of ‘5 to 1 hr, the temperature can 
to +0.02° C. Calibrated 
istantan thermocouples and a precision poten- 


stem down to to 


and to 


be 
iintained copper- 
meter are used for measuring the temperature 
the apparatus with an accuracy of +0.02° C. 
of the temperatures at the 
ts and outlets in the final saturator, A23, under 
ious conditions, indicate that at the inlets the 


Measurements 


air streams are on the average within 0.07° C 
me another, and at the outlets they are on the 


unidity Test Apparatus 


average within 0.02° C of one another. In the 
humidifying side of the final saturator, A23, the 
inlet and exit temperatures are on the average 
within 0.03° C of each other. The temperature 
in the test chamber is on the average 0.08° C 
warmer than the air emerging from the final 
saturator, A23. 

Since no further saturation can occur subsequent 
to the exit of the air from the final saturator, 
A23, the pressure at the exit may therefore be 
The difference 
between Ps and atmospheric pressure, P,, for 
various test chamber air velocities, is shown in 


the At 


considered as the Ps of eq 5. 


curve of air velocity 


figure 5. 


an 


15 
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Figure 5 Difference in pressure hetween the exit of the 
final saturator, P,, and atmospheric pressure, P4, for 
various test chamber air velocities 
of 525 ft/min, the assumption that P?.=)P, 


introduces an error in // of 0.4 percent, at 850 
ft/min, the error is 1.2 percent, and at 1,140 ft/min, 
the error becomes 2.5 percent. If eq 5 instead of 
eq 6 is used for calculating the relative humidity 
produced by the apparatus, these 


errors” are 


avoided. 
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IV. Calibration 
1. Theory 


In order to check the over-all performance of the 
apparatus, a series of tests was made to compare 
the actual humidity it produces with the humidity 
given by eq 6. The gravimetric method of mois- 
ture determination was selected as the means for 
measuring the water vapor content of the air in 
By 
test 


the test chamber under operating conditions. 


definition, the relative humidity in the 
chamber is 
Pu 


Ps 


H 


where py, is the partial pressure of the water vapor 
in the air, and p,, is the saturation pressure of the 
Assuming that air and water vapor 
60° C 
sure of 1 atmosphere obey Dalton’s law of partial 


water vapor 
from room temperature to and at a pres- 


pressures and the perfect gas law, then 


Pa+ Pu=Pe, 


ae 


VU. RT¢, 


and 
Mn 


paV Vl RT-, (10 
where p, is the partial pressure and 7, is the total 
pressure of the air in the test chamber, V is the 
volume, m, and m, are the weights of water and 
air respectively in volume V, MJ, and MV, are the 


molecular weights of water and air, F# is the uni- 


Borometer 


Gasmeter__ 


versal gas constant, and 7, is the absolut 
perature in the test chamber. 

Since the air in the air-water-vapor mixtu 
the test-chamber exhausts with negligible p 
drop into the atmosphere, then P-=P,, w! 
is atmospheric pressure. Dividing eq 9 by 
and substituting into eq 8 yields 

P, 

MoM » 

T im. M, 


Pw 


If, from a sample of air taken from thy 
chamber, the weight of water vapor and 
determined and the pressure measured, the 
pressure of the water vapor could easily by 
lated using eq 11, and the relative humidity y 
follow from eq 7. Figure 6 is a schemati: 
gram of the setup used for obtaining m 
P,. 


chamber and a sample of air withdrawn 


A glass probe was inserted into thy 


water vapor was removed from the sampk 


passing it through an absorption tube containing 


phosphorous pentoxide and accurately weighed§ 
The 


same sample of air were measured, and from thes@ 


volume, pressure, and temperature of 


data m, was obtained. 


Now 


Ihe 


Vd, 


where d, is the density of dry air at the pressur@ 


- 


and temperature of V,, the volume of dry aig 


In this particular setup, V, was measured by § 
vas meter at a 


laboratory wet temperature (g 


Woter 
Manometer 


4p Vacuum 


Valve 





Absorption 
Tube 
LI 


Mixing 


Soturated 
Air 


Ficure 6 


Ps 


FO” 


Borometer 


Vv 
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Diagram of gravimetric calibration setup. 
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smeter completely saturates the air passing 

it so that the partial pressure of the dry 
1/2Ap, 
’» Is the pressure at the inlet to the gas 


he gasmeter is given by P,—p, 
), is the saturation pressure of water vapor 
emperature, ¢,, of the air in the gas meter, 
is the pressure drop across the gas meter. 


ight of the dry air in the sample, in terms 


observable data and the density of dry air 


der standard conditions, becomes 


1/2A p)(0.001293) (273.1) 


— _ (13) 
(60(273.1- Sa) 


bstituting from eq 13 and 11 into eq 7 yields 
desired expression for the relative humidity in 
test chamber in terms of known and measur- 
e quantities 


P, P su 
V.(Ps—p, 2Ap)(0.00028898 ) 


] 
my(273.1 t t.) 


(14) 


1+ 


With the proportioning valve set to produce a 
lative humidity of 100 percent, and the pressure 
gulator adjusted to give an air velocity in the test 
amber of 525 ft/min, a number of gravimetric 
terminations was made in the temperature 
inge of 10° to 40° C. Then with the air 
ocity maintained at 525 ft/min, and at an 
bient temperature of —10°C, a run was made 


each setting of the proportioning valve. 
nally, at an ambient temperature of —10° C 
d with the proportioning valve set at 100 per 
nt, two runs were made, one at an air velocity 


850 ft/min and the other at 1,140 ft/min. 


2. Procedure and Results 


ln making a gravimetric determination, the 
midity apparatus was set into operation and 
ught to equilibrium conditions. The sample 
‘air was withdrawn by the glass probe at the 
te of 1 liter per minute. An attempt wes made 
colleet 50 mg of water vapor in the absorption 
«. This required withdrawing air for periods 
The weight of the water 
is determined with an accuracy of '; to 1 percent. 


irving from 's to 10 hr. 


For the duration of the run, frequent readings 
re made at equal time intervals of the tempera- 
res and pressures. The pressure, P,, at the in- 
tof the gas meter was determined with an ac- 

y of +0.2 mm Hg and the atmospheric pres- 
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sure, P,, was determined with an accuracy of 
t0.1 mm Hg, both by means of precision aneroid 
barometers. A mercury-in-glass thermometer, ac- 
curate to 0.1 deg, gave the temperature, f,, of the 
air in the gas meter. A water manometer that 
could easily be read to 's mm or better was used 
to obtain the pressure drop, Ap, across the gas 
meter. 

The results of these determinations are given in 
table 1, based upon the saturation pressures with 
respect to ice presented in the Smithsonian Tables," 
in Warmetabellen,’ and in Goff’s table.® It is esti- 
mated that the vapor pressures have been gravi- 
metrically determined with an accuracy of 1 to 2 
percent. The difference, D, without regard to 
sign, between each test result and the setting of 
The differ- 


ences depend upon the vapor pressure tables used 


the proportioning valve is also given. 
in the computation. The best agreement between 
the gravimetric determinations and the valve 
settings of the humidity apparatus is obtained by 
utilizing the Wirmetabellen values of vapor pres- 
TARLE 1 CGravimetric calibration 
Relative humidity 
satura " 
tion * * Smithso Wiirmeta 


temper nian bellen 
ature 


tesult D Kesult 


‘Smithsonian Meteorological Tables, Sth ed. (1939 
’ L. Holborn, K. Scheel, and F. Henning, Wirmetabellen der Physikalisct 
rechnischen Reichsanstalt (1919 
John A. Gotf, ASHVE Journal Section; Heating, piping, and air con 


ionin 155 (June 1945 
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sure. The average difference between the valve 
settings and the test results for the temperature 
range of —10° to —40° C and for air velocities 
up to 1,140 ft/min, based on the Warmetabellen 
is 1.7 pereent, on Goff’s table is 2.4 percent, and 
on the Smithsonian tables is 2.6 percent. These 
results show that the humidity apparatus, over 
the temperature and velocity ranges tested, pro- 
duced air having a relative humidity that aver- 
aged within 3 percent of the orifice ratio setting 
of the proportioning valve. 


V. Discussion 


This apparatus is particularly suited for the 
investigation and calibration of hygrometers, psy- 
chrometers, and dewpoint apparatus at tempera- 
tures below 0° C. The rotation of the propor- 
tioning valve and the adjustment of the pressure 
equalization valves are the only operations neces- 
sary for changing the relative humidity. The 
time involved is usually a matter of 1 or 2 seconds. 


This readily permits any device under test to be 


subjected to a sudden discontinuous change in 
relative humidity, a desirable advantage in the 
study of lag characteristics. The air velocity in 
the test chamber can easily be adjusted to any 
value up to 1,500 ft/min within several seconds 
by adjusting the pressure regulator. The size 
of an object that can be tested is limited by the 
inside diameter of the test chamber, which is 1'x 
in. For applications where high air velocities 


are not necessary, the diameter of the test ¢| 
may be increased to allow larger devices 
inserted. 

There is no provision in the apparatus fo 
ducing discrete and rapid changes in a: 
temperature. To change the ambient temp 
may require up to an hour to reach equi! 
conditions. 

There are two ways in which this apparaty 
be used to produce known dewpoints. If t! 
portioning valve is set to 100, then the en 
air will have a dewpoint equal to the test ch: 
temperature, or, more closely, to the tempr 
of the air in the final saturator. A chang: 
point is made by altering the bath temp 
On the other hand, if recourse to standard 
pressure tables and eq 4 is made, then for any g 
bath temperature, several vapor pressures, equiva 
lent to dewpoints at temperatures below that off 
the bath, are obtainable by rotation of the propor-] 
tioning valve. The dewpoints thus obtained may] 
have an average error of 0.3 deg C based or 
average error in the humidity of 3 percent 


The author thanks R. B. Kennard, for designing 
the proportioning valve and preliminary models.] 
and E. 


units, especially the saturators. 
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G. Clarke for constructing some ol 
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Gradual Damping of Solitary Waves 
By Garbis H. Keulegan 


This paper treats the problem of the damping by viscous action of translation waves. 


A shert exposition is given of Boussinesq’s boundary layer theory for wave motion, and 


expressions for the damping of rectangular and solitary waves are derived 


experimental results for solitary 


agreement is found to exist. 


Scott 


Russell's 


waves are compared with the theory, and satisfactory 


This fact makes it legitimate to apply the formulas developed 


to correct in model tests on harbors or in other tests of a like nature for the dissipative effects 


ir in shallow-water waves 


I. List of Symbols 


that occ 


wh,/H. 

breadth of channel. 

base of natural logarithms. 

rate of dissipation of energy of the wave. 

rate of decrease of energy of the wave 

energy in column of unit length and of width dz 

function of 

acceleration of gravity 

height of wave above undisturbed surface. 

initial height of wave above undisturbed surface 

depth of undisturbed water. 

constant defined by eq 35. 

constant defined by eq 44. 

length of wave 

kr. 

an integral defined by eq 41 

the pressure in the liquid. 

a variable. 

a variable, also the distance travelled by the 
wave. 

time, 

coordinate parallel to solid surface in direction 
of motion. 

coordinate parallel to solid surface at right angles 
to direction of motion. 

coordinate perpendicular to solid surface. 

velocity components of particles in the boundary 
layer parallel to z, y, z, respectively. 

velocity components of particles in the potential 
region parallel to xz, y, z, respectively. 

particle velocity in the potential region. 

z-component of body torce per unit mass due to 
gravity. 

J/-component of body force per unit mass due to 
gravity. 

z-component of body force per unit mass due to 
gravity 
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function of 
a variable, also a damping coefficient. 
a variable. 
a constant defined by eq 17. 
viscosity of the water. 
kinematic viscosity of the water. 
density of the water. 
a dissipation function. 
w=velocity of wave propagation. 


II. Introduction 


In the propagation of translation waves in still 
water, the main source of dissipation of the energy 
is the thin layer of liquid next to the solid bound- 
aries of the channel; that is, in the boundary layer 
The motion in the layer is laminar, and the de- 
termination of the velocities in it can be made by 
resorting to the ordinary equations of viscous 
motion. The motion of the liquid outside of the 
layer is potential, and there the values of the 
velocities are exactly those that would be obtained 
from the ordinary theory, assuming that the lami- 
nar layer is absent. The errors involved in the 
use of this method can be ignored, since the varia- 
tion of pressure in a cross section of the laminar 
laver is negligible. In fact, the gradual decrease of 
height of a solitary wave, the only wave of transla- 
tion, which under ideal conditions, travels without 
change of form, can be explained on the basis of 
this dissipation, at least, when the wave is of 


moderate proportions. The purpose of this paper 


is to develop a method of computing this energy 


dissipation and (to compare the results with ex- 
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perimental data as a means of verifying the 


accuracy of the theory. 


Ill. Dynamical Equation of the Laminar 
Layer 


The theory of the laminar layers with reference 
to wave motion was initiated by Boussinesq, and 
the essential elements of his method will now be 
given [1].! 

Consider the unsteady motion of a liquid near 
a flat solid boundary. For the sake of generality 
it will be supposed that the motion is two-dimen- 
sional. We select the rectangular system of axes 
r, y, 2, Where 2 is normal to the solid surface. The 
velocity components of particles along these axes 
At the solid bound- 
In the motion 


are u, v, and 0, respectively. 
ary, that is, at >=0, uw and ¢ vanish. 
of waves u, Ou/Or, and similar quantities are small, 
and hence their products may be neglected 
Owing to the thinness of the layers and their great 
lengths, only shears parallel to the solid boundaries 
are of consequence. Thus, the equation of motion 


for the liquid in the layer may be written as 


Ou lop, ‘ Oru 
ot por’ "dz? 
or l Op a -“) 
ot p oy ” O02? 
1 Op 
7 


abe Z, 


and 0 
where v is the kinematic VISCOSItY, p is the density 
of the liquid, p is the pressure, and Y, Y, Z are the 
components of the body forees per unit mass due 
to gravity. The third equation indicates that in 
the boundary layer the pressure p, the part due 
to disturbances, does not vary along the normals 
to the boundary. The pressure, p, in the viscous 
laver is the same function of z and y as in the 
potential region next to the viscous layer. 
Denoting the limiting values of u and ¢ at the 
boundary of the layer by uw» and 7, these last being 
the velocities in the potential region, we have for 
this region 
OUy 1 Op 
of por 
or 1 Op , 
dt ° by + Y. 


X,} 


Figures in bracket ndicate the literature references at the « 


paper 
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Subtracting these equations from thy 
sponding ones in eql, we obtain 


O(u— Uy O*(u— Uy) 


v 


and 


These are the basic boundary equations : 
rived by Boussinesq. 
they are applicable to a progressive wavy 
oscillation in a channel or to a standing wav: 
rectangular basin 

The limiting 
boundary conditions is introduced for conveniet 
Physically 1 


value, z=, appearing in 


in the mathematical analysis. 
boundary layer is very thin, and*the velocit 
u and ¢ attain the values wv and 7%, respective! 
finite distances away from the sol 
Moreover, if the waves examin 


at small 
boundaries. 
are progressive oscillatory waves moving in 
rectangular channel and the layer next to tk 


te 


bottom is considered, the axis of 2 is taken ver 
and the zs axis is in the same direction as | 
channel axis. For this case r=0, and w Is 
function of z, z, and ¢. When the boundary lay 
of one of the vertical walls is considered, the a 
of z is taken normal to the wall, the axis of » 
vertical, and, as before, the axis of x is parall 

the channel axis. Now wp, and 7, are the veloc 
in the potential region outside of the layer and a 
The case of | 


standing waves in a rectangular channel 0! 


functions of x, y, and ¢ only. 


rectangular basin is interpreted in_ th 
manner. 
As we shall be interested primarily u 


tion waves, we need consider only 


O(u— Uy) O*(u— Up) 
y ' , 


of O2= 
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In the form given abovg 


of 
ousl 
ripte 
ation 
trt7 


Wave 


the 
lining 
eent 
n this 


wing 


boundary conditions 


s the uniform velocity in the cross section 
function of x and ¢, assuming that the 


vels in the direction of x positive or z 


approximate theory of waves, if / is the 
vation of the wave surface, measured from 
sturbed water surface, and /7/ is the depth 
iter in the channel, for a wave traveling in 


tion of x positive: 


A 
-_— 


is the velocity of Wave propagation given 


ygll. (10) 


a wave of translation having any arbitrary 
but a small wave height, and for a solitary 
ndependently of height, for waves traveling 


he direction of « positive, we have 


wt 


wt). (12) 


he meaning of these expressions is that the 
e of waves does not change with translation. 
ously there are two methods available for the 
vription of the waves. We may observe the 
ition of the wave height with time at a fixed 
tr(r=0), or we may observe the variation of 
wave height with z at a fixed time t(f=0). 
Se select the second method of observation, and 

the main problem resolves itself into de- 
ining the velocities in the boundary layer 
e entire length of the wave at a given time, f. 
u this point of view, the basic equation now is, 
wing O/dt in eq 7 by —w(d/dr) 


Olu Us, vou Uy) 
Or O2 


the boundary conditions 
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In general, for ¢ constant, it will be assumed that 
Uy, Fi Bie { 15) 


IV. Distribution of the Velocities in the Vis- 
cous Layer. 


The expression 


é cos E a—dJ - ~ tb (16 
\ 2» 


where « and @ are independent of 2 and sv, as is 
well known, is a solution of eq 13 and satisfies the 
second equation of the boundary conditions, eq 14. 
Putting 2=0, and remembering that u=0 at 


0, we obtain 


r)). 17 


u COs [k a 


Thus we conclude that when the velocity of u 
in the potential region of the water is periodic and 
is given by eq 17, velocity « in the laminar layer 
is given by eq 16 

To obtain the velocities in the viscous layer 
when the wave consists of a long intumescence, 
we require a general solution with the value of 
given by eq 15. The general solution may be 
constructed from the elementary solution, eq 16 
Multiplying its right-hand member by F(a)da and 


integrating between the limits, + © and 2, We 


ae 
Fila ¢ COST Kla@ j " _ da, 
‘ | \ =| 


which still is a solution of eq 13 and satisfies the 


obtain 


second condition of eq 14. Since the permissible 


values of « are positive numbers, we finally take 


KW 
cos | Kla r) 2 da, 1S) 
\ Dp 


For, when we 


and this is the desired solution. 
make z=0, u=—0, we have 


F(a cos [x(a 


d \da, 


and, therefore, according to the rule of the Fourier 
F(r), 


first condition in eq 14 is also satisfied 


integral, it gives u indicating that the 
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Although eq 18 gives the distribution of veloci- 
ties in the boundary layer for the most general 
case, it is desirable to transform the solution into 
a form more directly manageable for calculation. 
That is, it is better to change the double integral 
into a single integral by effecting an integration 
with respect to x. With a view to making use of 


the two well-known definite integrals, 


: [" F (a)da [ é 


: F(a)da \° 


In the first double integral in eq 21 aSz, and this 


suggests the transformations 


and 


Introducing these transformations in the first 


double integral of eq 21, the integral changes to 


4 % F(a da . 


rT ¥ a a 


qs sds, 


In the 
second double integral of eq 21, aSz, and this 


and therefore vanishes by virtue of eq 19. 


suggests the transformation 


v q° 

and 

2w2* 

da 3 dq. 

vg 
Introducing these transformations in the second 
double integral of eq 21 and making use of eq 20, 
the final result ts 


) 


wr 


vq" 


ye | dq, 


| e-® cos (2s*+-qs)sds=0, 
e 
and 
. 


e~? cos (28° qs)sds 
) 


broken into two parts, as follow- 


ty! ). vd 


which is the desired solution and may be : 
with the Duhamel solution of the heat-« 
equation {2}. 


V. Loss of Energy in the Laminar Lay 


As it is assumed that the laminar lay: 
vicinity of the solid boundaries is thin and 
considerable length, it suffices to express 1! 
sipation function in the form 


ou 
od u ( a2 ), 
where yu is the viscosity of the liquid. 
rectangular column of indefinite height, the 
section being of unit length and of width ¢ 
noting the rate of dissipation of energy 
column by d/dt 6k), 


d i "2 (ou? 
dp OFA) pda | (=) dz 


* 


Now, since ue is independent of 2, 


("(= ; l [ O(u— Uy) |° 
dz 

Jo Oz Jo dz 

Integrating the right-hand member by 

have, since u— Up uy When z=0, and 


when z= ©, 


[- (2) dem m( 





u/Oz)o is the gradient of velocity at the 


i view of eq ia. 


32) d= wo (32),+5 ary 


ting in eq 23, 
bE) = udydruy ( ~ ) 4 
0 
pwd ydxr Da 


his expression is integrated between the 


hh 
ts +o and zs ©, the result gives the 
of energy dissipation in the entire length of 
laminar boundary layer. Denoting this loss 


i, dt, the integration gives 


te siyfw(B) as, ev 


at the infinite distances, z © and « 2, 
velocities uw and wp both vanish. This formula 
om Boussinesgq [1]. The interpretation of the 
ula is that for the determination of the rate 
ssipation of energy in the entire wave it is 
cient to consider the gradient of the laminar 

ties at the solid boundary and the velocities 
he region outside of the boundary layers. 
his result can be extended to limited waves 

a new interpretation. Let the posterior and 
rior feet of the wave be at the points x, and 2x», 
ectively. Since the thickness of the viscous 
rat the front of the wave; that is, at point 2», 

gligibly small, the expression for the rate of 
sipation of energy for the limited wave per 

width of channel may now be written, since 


. ndependent of 2, as 


0 Se) dru [ae f° (3" 


“, is the laminar velocity at the point 2. 
interpretation of the two integrals on the 
-hand side of the equation is as follows. The 
ntegral (double) represents the rate of energy 
‘pation in the laminar boundary layer under 
wave. The second integral represents the 
tic energy left behind in the channel as the 
travels onward. As far as the wave is 
erned, the energy left behind is lost, so that 
n the case of a limited wave the integral on 
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the left-hand side represents the rate of loss of 
energy of the wave per unit width of channel. 


VI. Gradual Extinction of Rectangular 
Waves 


The rectangular wave constitutes the simplest 
example of a translation wave, and the application 
of the above theories to obtain the expression for 
the gradual damping of the wave with distance 
involves no difficulty. In a rectangular wave the 
superelevation is constant, and the wave length 
is finite. 

Denote the constant superelevation of the wave 
by Ay, the length by L, and the particle velocities 
in the potential region by V. Assuming that the 
wave is traveling in the direction of positive z, 
we have 

w 7 w= 9H. (25) 
Suppose that at time f=0, the front of the wave 


is at the station c=Z (see fig. 1). Then for the 











FIGURE 1. Surface profile of rectangular wave 


description of the function F appearing in eq 15, 
we have 
Fiz 0), 


F(z)=V, 
F(x)=0, 


») 


Accordingly, using eq 22, we find the boundary 


laver velocities to be 


e~* dB, 


and 


sy rr 
e~*dB+ 
VyUdJz 


~ / o. 
2V =a 
when Since, in the present example, the 
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rate of dissipation of the energy of the wave per 


unit width of channel reduces to 


dk, “d Ou - 
= a wo(,.) dz, (29) 
by the use of eq 23, it is sufficient to consider only 
eq 27. Differentiating the latter equation with 
respect to z and remembering that V is independ- 
ent of 
ou V i "= 
O- y rv L—sr)/w 
Putting ° 0, the velocity vradient at the solid 


boundary is 
Ou 
\d2).= VaniL 


Substituting in eq 24 and remembering that 
u V, 

dk: (" Vder 

at \ 7)  f J Ww 
Integrating and afterwards substituting from eq 25, 


we have finally 


dk 2 


dt pd WH h 4%, “yp. oO 


\7 


which is the value of the energy lost per unit time 


and per unit width of the channel. As was ex- 
plained previously, the loss is made up of two 
One 


energy in the boundary layer under the wave; the 


parts part represents the dissipation of 
other part represents the kinetic energy left in the 
liquid and in the vicinity of the solid boundary 


as the wave leaves a given station. 


The energy of the wave, half potential and half 


kinetic, is given by 


E,—gph?L, 


when computed for unit width of channel. The 
variation of the energy with the time is 
dE. dh, 
“==2gph,L ,- (31) 
dt 29 ph yt at | 


Denoting the distance traveled by the wave by s, 


dh, dh, ds 
dt ds dt’ 
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and as ds/dt represents the velocity of the 


dh, dh 
dt vgH 


ds 
Substituting this in eq 31, 


dh 
dt d 


18 


2pg *H] 7 Lh 


By the principle of the conservation of 
the decrease of the wave energy per u 


must equal the dissipation of energy | ; 


time. Accordingly the right-hand mem}» 
30 and 32 are equal to each other in absolute va 


With the 


of the terms after the comparison is nh 


but they are of opposite sign 


have 
/ / l/l 
(77) “y) 07, V gene yt 
which is the differential equation giving e 
of decrease of the wave height with the dist ® 


of travel. The equation was derived by assu 
that the channel is very wide and thi 


If the ch; 


is narrow, however, the second term in the 


the side walls could be neglected. 
hand member must be multiplied by (1 


where B denotes the width.of the channel. 1 
is, in a channel of finite width, 


h,\", (hy 
(Ci) “Gi ~ 74 


=a 7) (Gen .) (41) 


The solution of the equation is, subject , 
assumption, h=ho, when s=0, 
h, ho “UW % 
é 
ll H i 


where 
| 2H\ (L\? y 4 

. ve i4 B a7) (, spp) ' 

A 


Accordingly, the decrease of wave he 
distance is given by an exponential, and a 
rithmic decrement that is independent of the hei 
of the wave exists. 
of constant depth, the damping decreases 


For waves traveling in wag 
increase of wave length, with the increas 
channel width, and with the decrease of k 


Viscosity. 
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Vl. iradual Extinction of Solitary Waves 


Denoting the superelevation of the crest by hy, 
superelevation of the wave profile for f con- 
(0) is given approximately [3] by 


h=h, sech? ke, 


1 /3hy (26 
36) 

HV 4H : 

dimensionless plot of the wave profile, when 
depth of water is twice the crest elevation, 


shown in figure 2. Using eq 9 and 12, we may 





Surface profile of solitary wave, H/2). 
te for the particle velocities in the potential 
sion the expressions, assuming that the wave is 
veling in the direction of s negative, 


A sech? kr; A=w ri 


he region of the viscous boundary layers, using 


®dgB (38) 


2A kw2 
- ; = |, sech* ( kat 43? ) é 


ferentiating with respect to z and remembering 
fy 1S independent of 2, 


Ou 4A 
0 \ 7. 


( sech?( ket cca ) tanh ( kat 
ay kw2 Pde. 


4vp? } 2vp? ‘ 


limiting value of the integral when 2z ap- 
aches zero cannot be obtained directly, and a 


sformation of variables becomes necessary. 


and 
kw 2 


V 4y B 


a 
so that 
\ * A dp da, 


we have 


kw . » 
seclh? (hi 
0 


Ou 44 
Oz va an 


+e?*) tanh (n- 


, 4 < 
a’ je " da. 
This gives for z=0 


4A kw ‘ > 
sech? (n-+-a’*) tanh (n 
\7 \ v Jo 


tar) da. 


(39) 


_ 
rhe 


integral appearing in the right-hand side is diffi- 


rigorous determination of the definite 


cult, and sufficiently accurate results can be 
obtained if the indicated integration is replaced 


by the process of summation. Writing 


V(n) sech? (n+ a’) tanh (n+-a*)da, (40) 


the desired summation for a specific (nm) was car- 
ried out by assigning to the interval da the value 
0.1. The results of such determinations are given 
in table 1. 
Since ¥(n) is proportional to the velocity gradient 


The graph of }(n) is shown in figure 3. 


at the channel bed, the curve may be interpreted 
as giving the variation of this velocity gradient 
The velocity gradient at the 


along the wave. 
The maxi- 


O.35. 


wall reverses its sign at n 


mum value of the velocity gradient at n 


TABLE 1 Vumeri-al values of the function Yin 














n 


Fiacure 3 Variation of velocity gradve nt at wall with 


lance from center of wave. 


For the case when the wave height is half of the 
undisturbed water depth, the reversal occurs at 

0.9/7. The maximum occurs at 
0.57/17; that is, 
position is given by 2=U of 


the point, 4 
the point, 7 ahead of the wave 
crest. The crest 
R=. 

In view of the substitutions from the eq 37, 39, 
and 40, and the transformation, kr=n, the formula 
for the rate of dissipation of energy, eq 23, reduces 
to 


im Nt) sech? ndn, 


dE, 4A? ke | 
dt yrk\ Zz 


which may be written also as 


dk, 4A? WY ,, 
dt \7 PX k : 


where 


(- Yn 


The method of summations, 


sech? ndn. 


selecting n=0.1, 
gave for N the value 0.316, or approximately 1/7. 
Using the value of & from eq. 36 and A from eq. 37 
ygll, the above ex- 
pression for the rate of dissipation of energy per 


and remembering that 


unit width of channel becomes 


dk 4 ly4 

yl /2 5/474 (49) 

dt Sr (s) pv Pepa, 42 

As regards the energy of the solitary wave for 
unit width, we have 

So . 

" g( Hh, 8” 

) 


By: 





Differentiating with respect to ft, and remen, be) 
that 


a arth 


we find for the rate of decrease of the e: 
the wave 

dk, 4 . dh 

, eg? IPhY? |! 
dt y3 ds 
By the principle of conservation of energy 

right-hand members of the eq. 42 and 43 are equ; 
in absolute value, but they are of opposite sig 
Making the comparison and rearranging the t 
we obt iin 


h, : h, l v , 
Ci) an) =~ 3y gee (y) 


This equation is derived by assuming that 

channel is very wide, and, therefore, the effect | 
the vertical side walls can be neglected 
however, the channel is narrow, the second | 
in the left-hand member must be multiplic: 
{1+ (211 B)), where B denotes the width of ¢! 


nel. That is, in a channel of finite width 


v 


h, . h, ] 2T1\ Q 
(7) d(7,)+3('+ BV gemey) 


The integration gives, denoting the initial 
of the crest height by Ajo 


h, hc\'* —e 
(i (77) —Kip 


where the coefficient of damping has the val 


(42H! 


A=i9 BN g@t 


foll 
This gives the variation of the crest height nlol 
solitary wave with the distance traversed eet 
Comparing eq 33 and 44, it is seen tha 
tinction of a limited wave takes place in a 
very different than it does with an unl 
wave. Although in the former case th 
tion is characterized by a logarithmic dl 
independent of height, this is not tru 
solitary wave. The reason for this is to 
in the fact that the effective wavelength 


tary wave is a function of the crest heigh 
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rease of the crest height, the effective 
nereases [4]. We may now accept the 
statement that a logarithmic decrement 
dent of the wave height exists only for a 
at Maintains a constant wavelength during 
on. Thus the damped oscillations of a 
na U-tube, of a standing wave in a rec- 
r basin and the damping of progressive os- 

wave are characterized by a logarithmic 
ent independent of the magnitude of the 


rons. 


John Scott Russell's Experiments on 
Solitary Waves 


test’ the validity of the theoretical results 
ve given, John Scott Russell’s experimental 
ta on solitary waves will be considered [5]. Al- 
ugh his tests were made to determine the law of 
ocity of propagation of solitary waves, the 
resentation of data is such that they can be 
wd also to determine the law of damping of 
ves. The experiments were made more than 
entury ago in a wooden rectangular channel, 


wide and 20 ft. long. Altogether some 55 
sts are shown for various water depths and for 
ous methods of wave generation. 
Regarding the tests, Scott Russell states: ‘‘ The 
yinal experiments are themselves given at the 
of this paper for the purpose of enabling 
rone who may be disposed to make use of them 
any future purpose, either of framing or testing 
heory, to make use of them much in the same 
vas if he had himself made the experiments.” 
farsightedness and disinterestedness of the 
stigator is evidence of a scientific attitude of 
highest type and are cited by the present 
ter, to whom Russell's data have proven in- 
ble, as a model that other investigators might 
follow. 
Cnfortunately the temperature of the water of 
tests Is not stated Because of this, the de- 
(d comparison between theory and experiment 
be made by deducing the temperatures of 
ter from the observed damping of the waves 
comparing these values with the prevailing 
temperatures on the dates of the experiments. 
in consulting the original data, it was seen that 
tests could be grouped together and average 
wes obtained for individual water depths /7. 
ethod of arriving at average values may be 


adual Damping of Solitary Waves 


explained by considering the measurement for //= 
17.9 em. It was thought that five individual 
tests with this depth ought to be sufficient. The 
original data furnish values of /', the height of the 
wave, and s, the distance travelled. From these, 
the ratios A'//7 and s/H/ were formed and plotted. 
Since in these tests the initial heights hj were not 
the same, the plotting gave separate curves. 
Keeping one of the curves fixed, the others were 
shifted by different amounts parallel to the axis of 
abscissas so as to give a more or less congruent 
pict ure. The positions of the observed points after 
the shiftings are shown in figure 4. The centers of 
gravity of the points within selected increments of 
the abscissas were obtained and are given in table 
2. Similar determinations were made for other 


depths and these are given likewise in table 2. 


TABLE 2 Eatinction of solitary wares; average values from 


R tl’ 
TuUSse s experiments 


The data for /7=17.9 em from table 2 are 
plotted in figure 5. The smooth curve drawn 
through the points gives the adjusted values of 
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FiGureE 4 
the wave heights. The quantities A,///7 and s/// 
Values of the 


and 


were read from the smooth curve. 
fourth Hh determined 
plotted in figure 6 with a convenient shift in 
the origin of s/7/ 
2 were treated in a similar manner, and the re- 


root of were are 


The remaining data from table 


sulting curves are shown in the figures 6 and 7. 
The distribution of the points is linear, a cir- 
The 


slopes of the lines determine the coefficients of 


cumstance that accords very well with eq 44 
damping, K The viscosities of the water can be 
obtained from eq 45, using the observed values of 
K, the corresponding values of //, and putting B, 
the width of flume, equal to 30.48 em. The indi- 
vidual determinations of water temperature are 
shown in table 3, and the average value is found 
to be 15.8° C 
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Data from Scott Russell's experiments on solitary waves 
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Figure 5. Decrease in height of a solitary wave 
traversed 


H=17.9em 
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FicuRE 6. Coefficient of damping of solitary waves 


H1=2.69 cm, K=9.0x10 . 1=7.61 em, K=4.85x 104 


. 1=17.90 em, K=3.85x10 








1000 1500 2000 2500 3000 3500 


s/t 
('o¢ fhic rent of damping of solitary waves 


1.67x10- . 1=12.92 em, K=4.25.10' 


I1=15.80 em, A 


10.02 em, A 


4.25x 104 


Determination of vater temperatures in Russell's 


er periments 


0. 01168 
O1046 
O1133 
O1132 
O1215 


O1048 


This average temperature of 15.8° C can be 
mpared with the prevailing mean air tempera- 
res at the time of the experiments. The Russell 
periments were made in kdinburgh during the 
riod, August 3 to 8, 1837. It may be said that 
en water has been introduced from an outside 
vervoir into a wooden flume of small depth in- 
¢ a laboratory, during the latter months of 
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summer, the working temperature of the water 
of the flume is best approximated by the morning 
air temperature, rather than the minimum or 
maximum daily air temperatures. 
courtesy of the Meteorological Office, Air Ministry, 
Edinburgh, the air temperatures for the dates 
corresponding to Russell’s experiments were com- 


Through the 


municated to the Liaison Officer, Meteorological 
Office, U. S. Weather Bureau, Washington, D. C. 
The data later were transmitted to the author by 
the Chief of the Division of Climatological and 
".S. Weather Bureau. The 
author wishes to express his appreciation of these 
favors. The Royal Society of Edinburgh have 
permitted a copy to be made of the pertinent part 
of Mr. Adie’s original register for 1837, and the 
part of the data relating to the temperatures is 
It is the opinion of the Edin- 


Hydrologic Services, | 


given in table 4. 
burgh Meteorological Office that the morning and 
the evening readings of the dry bulb were taken 
at 10a.m.and10p.m. Accordingly, the temper- 
atures described as the morning temperatures 
should be close to the water temperatures of the 
The average of the morning 
temperature for the 7 days is about 14.6°C. This 
agrees reasonably well with the average temper- 
atures obtained from the analysis of the wave 
ia CC. 


Russell experiments. 


data; that is, 
Edinhurgh temperatures, Auqust 1837 


TARLE 4 


Mir Max 


Morning Evening 
mum mum 


The favorable agreement between the two tem- 


peratures and the fact that the inverse of the fourth 


root of the wave height varies linearly with the 
distance of the travel of the waves are proofs 
enough to show that the present theory of viscous 
action is sufficient to evaluate the wave dissipa- 
The next 
step in the investigations of this type should be 


tion for the water depths considered 


to determine what depths of water will result in a 
type of dissipation other than that caused by 


Viscosity. 
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IX. Application to Model Harbor Tests 


To the extent that the theoretical formula for 
the damping of solitary waves is satisfactory, at 
least when the waves are moving in shallow water, 
we can apply the results of the above theory to 
model tests of wave action. The models are of 
the type in which the depths of water involved 
are small and in which the pattern of the individual 
waves can be approximated as a succession of 
solitary waves. No serious error is made by re- 
placing the ordinary pattern of the shallow-water 
oscillatory waves by a succession of solitary waves 
if this is done for the purpose of evaluating the 
energy dissipation. The method has the advan- 
tage that the wavelengths do not enter into the 
computations, so that the comparisons may be 
made on the basis of wave height only. 

One problem that most readily suggests itself 
and is susceptible to approximate analysis is the 
model test of wave action ina harbor. Ordinarily 
the heights of the waves approaching a model 
harbor are determined on the basis of Froudian 
scales after the respective heights of the prototype 
waves have been selected. For the actual model 
these values of the initial wave heights must be 
increased over those determined by the Froude 
law in order to correct for the energy dissipation 
in the model harbor. In the prototype the dissi- 
pation is not large, so that it can be ignored, and 
this fact facilitates the determination of the proper 
model wave heights The exact determination of 
the desired model wave heights in the approaches 
to the model harbor depends on the length and 
depth proportions of the harbor. The evaluation 


of the energy dissipation may be made sjep }y 
step for various consecutive segments of the moc, 


The average of the time rate of dissipation {or | 
entire harbor area can then be evaluated 7) 


’ 


height of the approaching Waves is increased 


such an extent that the increased flow of 


just balances the energy dissipation in the harbor, 


The author expresses his appreciation to Herbe; 
N. Eaton, Chief of the Hydraulics Section, for | 
editorial review of the paper and for constructiy 
suggestions; also to George W. Patterson, | 
former colleague, for advice and the review of | 
mathematical methods used in the paper 
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